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One of the most striking results from the Relativistic Heavy Ion Collider is the strong 
elliptic flow. This review summarizes what is observed and how these results are com- 
bined with reasonable theoretical assumptions to estimate the shear viscosity of QCD 
near the phase transition. A data comparison with viscous hydrodynamics and kinetic 
theory calculations indicates that the shear viscosity to entropy ratio is surprisingly 
small, r)/s < 0.4. The preferred range is rj/s ~ (1 ^ 3) X l/47r. 



1. Introduction 

One of the most striking observa tions from the Relativistic Heavy Ion Collider 
(RHIC) is the very large elliptic flowffl^. The primary goal of this report is to explain 
as succinctly as possible precisely what is observed and how th e shear viscosity can 
be estimated from these observations. The resulting estimates l -J l ^P I OM 
indicate that the shear visc osity to entropy ratio jy/s is close to the limits suggested 
by the uncertainty princ iple^, and the result of A/" = 4 Super Yang Mills (SYM) 
theory at strong couplin g I I 

s An 

These estimates imply that the heavy ion experiments are probing quantum kinetic 
processes in this theoretically interesting, but poorly understood regime. Clearly 
a complete understanding of nucleus- nucleus collisions at high energies is extraor- 
dinarily difficult. We will attempt to explain the theoretical basis for these recent 
claims and the uncertainties in the estimated values of rj/s. Additionally, since 
the result has raised considerable interest outside of the heavy ion community, this 
review will try to make the analysis accessible to a fairly broad theoretical audience. 



1.1. Experimental Overview 

In high energy nucleus-nucleus collisions at RHIC approximately 7000 par- 
ticles are produced in a single gold-gold event with collision energy, ^/s — 
200 GeV/nucleon. Each nucleus has 197 nucleons and the two nuclei are initially 
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Fig. 1. Overview of a heavy ion event. In the left figure the two nuclei collide along the beam axis 
usually labeled as Z. At RHIC the nuclei are length contracted by a factor of 7 ~ 100. The right 
figure shows the collision vertex of a typical event as viewed in a schematic particle detector and 
shows a few of the thousands of charged particle tracks recorded per event. The angle 6 is usually 
reported in pseudo-rapidity variables as discussed in the text. 



length contracted by a factor of a hundred. The transverse size of the nucleus is 
Rau ~ 5 fm and the duration of the event is roughly of order ~ Rau/c. Fig. [T]shows 
the pre-collision geometry. Also shown is a schematic of the collision vertex and a 
schematic particle detector. 

Usually the two nuclei collide off-center at impact parameter b and oriented at 
an angle '^rp with respect to the lab axes as shown in Fig. [2] During the collision 
the spectator nucleons (see Fig. [2]) continue down the beam pipe, leaving behind 
an excited almond shaped region. The impact parameter b is a transverse vector 
b = {bx, by) pointing from the center of one nucleus to the center of the other. As 
discussed in Section [5] both the magnitude and direction of b can be determined on 
an event by event basis. We will generally work with reaction plane coordinates X 
and Y rather than lab coordinates. 

The elliptic flow is defined as the anisotropy of particle production with respect 
to the reaction plane (see Fig. [5] and Fig. [3]) 

\Px+Py/ 

or the second Fourier coefficient of the azimuthal distribution, (cos(2(0 — ^rp))). 
Elliptic flow can also be measured as a function of transverse momentum px = 
Y^p^ -|- py by expanding the differential yield of particles in a Fourier series 

I dN 1 rfA^ , „ , ^ , , ^ 

— , , ——{1 + 2v2{pt)cos2{<P-^rp) + ...) . (2) 

Pt dyapTdq) ZirpT dyapT 

Here ellipses denote still higher harmonics, V4, vq and so on. In addition the flow 
can be measured as a function of impact parameter, particle type, and rapidity. For 
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Fig. 2. A schematic of the transverse plane in a heavy ion event. Both the magnitude and direction 
of the impact parameter b can be determined on an event by event basis. X and Y label the 
reaction plane axes and the dotted lines indicate the lab axis. 'I'^p is known as the reaction plane 
angle. 



a mid-peripheral collision {b ~ 7fm) the average elliptic flow {V2) is approximately 
7%. This is surprising large. For instance, the ratio of particles in the X direction 
to the y is 1 + 2v2 ■ 1 — 2v2 — 1-3 : 1. At higher transverse momentum the elliptic 
flow grows and at px ^ 1-5 GeV elliptic flow can be as large as 15%. 



1.2. An Interpretation of Elliptic Flow 

The generally accepted explanation for the observed flow is illustrated in Fig. [3l 
Since the pressure gradient in the X direction is larger than in the Y direction, the 
nuclear medium expands preferentially along the short axis of the ellipse. Elliptic 
flow is such a useful observable because it is a rather direct probe of the response 
of the QCD medium to the high energy density created during the event. If the 
mean free path is large compared to the size of the interaction region, then the 
produced particles will not respond to the initial geometry. On the other hand, if 
the transverse size of the nucleus is large compared to the interaction length scales 
involved, hydrodynamics is the appropriate theoretical framework to calculate the 
response of the medium to the geometry. In a pioneering paper by OUitrualt, the 
elliptic flow observable was proposed and analyzed based partly o n the c onviction 
that ideal hydrodynamic models would vastly over-predict the 

However, calculations based on ideal hy drodynarnics do a fair to reasonable job 
job in rep roduci ng the observed elliptic flo\\ l'*^^ l ^^ l ^^ l ^^ l ^^ . This has been reviewed 
elsewhere l ^3 | 24 | Nevertheless, the hydrodynamic interpretation requires that the 
relevant mean free paths and relaxation times be small compared to the nuclear 
sizes and expansion rates. This review will assess the consistency of the hydro- 
dynamic interpretation by categorizing viscous corrections. The principle tool is 
viscous hydrodynamics which needs to be extended into the relativistic domain 
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Fig. 3. The conventional explanation for the observed elliptic flow. The spectators continue down 
the beam pipe leaving behind an excited oval shape which expands preferentially along the short 
axis of the ellipse. The finally momentum asymmetry in the particle distribution V2 reflects the 
response of the excited medium to this geometry. The dot with transverse coordinate x = [x, y) 
is illustrated to explain a technical point in Section |2] 



in order to address the problems associated with nuclear collisions. This problem 
has received co nsiderable atte ntion recently and progress has b een achieved both 
at a conceptual I^5l26l27l28l29]g^j^^ practical level ISEIlQlTOlSTIIl] Generally macro- 
scopic approaches, such as viscous hydrodynamics, and microscopic approaches, 
such as kinetic theor y, are converging on the implications of the measured elliptic 
flow | 29 | 32 | 12 | 33 | 6 | 34 | 'pjjgj-g has never been an even remotely successful model of 
the flow with 77/s > 0.4. Since rj/s is a measure of the relaxation time relative to 
h/ksT (see Section [S]), this estimate of 77/s places the kinetic processes measured 
at RHIC in an interesting and fully quantum regime. 



2. Elliptic Flow Measurements and Definitions 

The goal of this section is to review the progress that has been achieved in mea- 
suring the elliptic flow. This progress has produced an increasingly self-consistent 
hydrodynamic interpretation of the observed elliptic flow results. This section will 
also collect the various definitions which are needed to categorize the response of 
the excited medium to the initial geometry. 



2.1. Measurements and Definitions 

As discussed in the introduction (see Fig. [2|) both the magnitude and direction of 
the impact parameter can be determined on an event by event basis. The magni- 
tude of the impact parameter can be determined by selecting events with definite 
multiplicity for example. For instance, on average the top 10% of events with the 
highest multiplicity correspond to the 10% of events with the smallest impact pa- 
rameter. Since the cross section is almost purely geometrical in this energy range 
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this top 10% of events may be found by a purely geometrical argument. This line 
of reasoning gives that the top 10% of events are produced by collisions with an 
impact parameter in the range 

nb ^ 

0<b<K, where 10% = — ^ , (3) 

Ctot 

and (Ttot — 7r(2i?yi)^ is the total inelastic cross section. After categorizing the top 
10% of events we can categorize the top 10-20% of events and so on. The general 
relation is 

/ , \ 2 

2— j ~% Centrality. (4) 

Here we have neglected fluctuations and many other effects. For instance there is a 
very small probability that an event with impact parameter 6 = 4 fm will produce 
the same multiplicity as an eve nt with b = fm. A full discussion of these and many 
other issues is given in The end result is that the magnitude of the impact 

parameter b can be determined to within half a femptometer or so . 

Now that the impact parameter is quantified, a useful definition is the number 
of participating nucleons (also called "wounded" nucleons). The number of nucleons 
per unit volume in the rest frame of the nucleus is pa{^ — Xq, z), were x — Xq is the 
transverse displacement from a nucleus centered at Xq , and z is the longitudinal di- 
rection. These distribut ions are known experimentally and are reasonably modeled 
by a Woods-Saxon forrrPS. The number of nucleons per unit transverse area is 

/oo 
dzp^(x-Xo,z). (5) 
-C30 

Then, after reexamining Fig. [3l we find that the probability that a nucleon at 
X = (a;, y) will suffer an inelastic interaction passing through the right nucleus 
centered b/2 ^ (+b/2, 0) is 

l-exp(-CT™rA(x-b/2)) , 

where (Tnn — 40 mb is the inelastic nucleon-nuclcon cross section. The number of 
nucleons which suffer an inelastic collision per unit area is then 

^ - r4(xi + b/2) [1 - exp (-(7«,T^(xx - b/2))] 



dxdy 

+ rA(xi-b/2)[l-exp(-a™T^(xx+b/2))] . (6) 

Finally, the total number of participants {i.e. the the number of nucleons which 
collide) is 

For a central collision of two gold nuclei the number of participants Np ~ 340 nearly 
equals the total number nucleons in the two nuclei, N = 394, leaving about fifty 
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Fig. 4. The standard Glauber eccentricity es,part as a function of the number of participants. 
jYmax ^ g^Q ^j^g maximum number of participants in a central AuAu event and Ra — 6.3 fm is 
the gold radius. The top axis shows the translation between impact parameter and participants. 
The root mean square radius -Rrms and the standard Glauber eccentricity are given in Eq. 1)13^ 
and Eq. JTSt . 



spectators. By comparing the top axis in Fig. [4]to the bottom axis, the relationship 
between participants, impact parameter 6, and centrahty can be determined. 

The reaction plane angle, rp, is also determined experimentally. Here we will 
describe the Event Plane method which is conceptually the simplest. Assume first 
that the reaction plane angle is known. Then the particle distribution can be ex- 
panded in harmonics about the reaction plane 

dN 

— cx l + 2i;2cos(2(0-*i^p)) + ... (8) 
acp 

If the number of particles is very large one could simply make a histogram of the 
angular distribution of particles in an event with respect to the lab axis. Then 
the reaction plane angle could be determined by finding where the histogram is 
maximum. This is the basis of the event plane method'^. For all the particles in 
the event we form the vector 

Q={Qx, Qy ) - COS 20, , ^ sin 2cl)}j . (9) 
Using the continuum approximation, Qx — J d(f> dN / d(f> cos{2(f>) , we can estimate 
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the reaction plane angle "if bp, from the Q-vector 

= (cos(2*2), sin(2«'2)) ~ (cos(2*flp) , sin(2*flp)) . (10) 

1^1 

Then we can estimate the elliptic flow as ^2'''^ ~ (cos(2(0i — ^'2)))- The estimated 
angle 5*2 differs from '^rp due to statistical fluctuations. Consequently vf^^ will 
be systematically smaller than V2 since ^'2 is not ^'pp. This leads to a correction 
to the estimate given above which is known as the reaction plane resolution. The 
final result, after considering the dispersion of "if 2 relative to the true reaction plane 
angle 5'pp is 

^obs 

W2 = ^ where ^ = (cos2(*2 - *flp)> ■ (11) 

In practice the resolution parameter 1% is estimated by dividing a given event into 
sub-events and looking at the dispersion in '^2 between different sub-events. 

There is a lot more to the determination of the event plane in practice. Fortu- 
nately the various methods have been reviewed recently^ZI. An important criterion 
for the validity of these methods is that the magnitude of elliptic flow be large 
compared to statistical fluctuations 

4 » ^ • (12) 

For V2 — 7% and N ~ 500 we have ~ 2.5. Since this number is not particularly 
large the simple method described above is not completely adequate in practice. The 
resolution parameter is ^ ~ 0.7 in the STAR experiment. At the LHC, estimates 
suggest that the resolution parameter could be as large as'^^ ^ ~ 0.95. Current 
methods use two particle, four particle, and higher cummulants to remove the effects 
of correlations and fluctuations. These advances are discussed more completely in 
Section [2.31 and have played an important role in the current estimates of the shear 
viscosity. The current measurements provide a unique theoretical opportunity to 
study systematically how hydrodynamics begins to develop in mesoscopic systems. 

We would like to measure the response of nuclei to the geometry. To this end, 
we categorize the overlap region with an asymmetry parameter e^.part 



ts.part — ^^2 



(13) 



Traditionally the average (. . .) is taken with respect to the number of participants 
in the transverse plane, for example 

(,^^.»> = ±/d.d,(,^-.^)il^. (14, 

We will explain the "s,part" label shortly; for the moment we return to Fig. IH 
which plots the asymmetry parameter versus centrality and also shows the the root 
mean square radius 
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which is important for categorizing the size of viscous corrections. 
2.2. Interpretation 

We have collected the essential definitions of e, centrality, and V2, and are now 
in a position to return to the physics. The scaled elliptic flow W2/e measures the 
response of the medium to the initial geometry. Fig.[5]shows V2{pt)/£ as a function 
of centrality, 0-5% being the most central and 60-70% being the most peripheral. 
Examining this flgure we see a gradual transition from a weak to a strong dynamic 
response with growing system size. The interpretation adopted in this review is 
that this change is a consequence of a system transitioning from a kinetic to a 
hydrodynamic regime. 

There are sever al theoretical curves based upon calculations of ideal 
hydrodynamic j ^^ l ^^ l which for px < 1 GeV approximately reproduce the observed 
elliptic flow in the most central collisions. Since ideal hydrodynamics is scale invari- 
ant (for a scale invariant equation of state) the expectation is that the response V2/£ 
of this theory should be independent of system size or centrality. This reasoning is 
borne out by the more elaborate hydrodynamic calculations shown in the figure. 
On the other hand, the data show a gradual transition as a function of increasing 
centrality, rising towards the ideal hydrodynamic calculations in a sy stematic way. 
These trends are captured by models with a finite mean free patlt^. 

The data show other trends as a function of centrality. In more central collisions 
the linearly rising trend, which resembles the ideal hydrodynamic calculations, ex- 
tends to larger and larger transverse momentum. We will see in Section [5] that 
viscous corrections to ideal hydrodynamics grow as 




where L is a characteristic length scale. Thus these viscous corrections restrict the 
applicable momentum range in hydrodynamic^. In more central collisions, where 
£„fp/L is smaller, the transverse momentum range described by hydrodynamics 
extends to increasingly large px- These qualitative trends are reproduced by the 
more involved viscous calculations discussed in Section [HI 

To conclude this section, we turn to Fig. [6] which compares the elliptic protons 
and pions to the flow of the multi-strange hadrons f2~ and 4>. (These hadrons have 
valence quark content sss and ss respectively.) The important point is that the n~ 
is nearly twice as heavy as the proton and more importantly, does not have a strong 
resonant interaction analogous to the A. For these reasons the hadronic relaxation 
time of the il~ is expected to be much longer than the duration of the heavy ion 
evenl!32l. Nevertheless the fl shows nearly the same elliptic flow as the protons. This 
provides fairly convincing evidence that the majority of the elliptic flow develops 
during a deconfined phase which hadronizes to produce a flowing D,~ baryon. 
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Fig. 5. Elliptic flow V2{pt) as measured by the STAR coUaboratioJSSHOl for different central- 
ities. The measured elliptic flow has been divided by the eccentricity — ehydro = es,part in this 
work. The curves are ideal hydrodynamic calculations based on Refs.EMU rather than the viscous 
hydrodynamics discussed in much of this review. 



2.3. The Eccentricity and Fluctuations 

Clearly much of the interpretation of elliptic flow relies on a solid understanding of 
the eccentricity. There are several issues here. First there is the theoretical uncer- 
tainty in this average quantity. For example, so far we have defined the "standard 
Glauber participant eccentricity" in Eq. (fT3|) . An equally good definition is pro- 
vided by collision scaling. For instance, one measure used in heavy ion collisions is 
the number of binary nucleon-nucleon collisions per transverse area 

^J^ = a™r^(x-Kb/2)T^(x-b/2), (17) 
axay 

Then the eccentricity is defined with this A'coii weight in analogy with Eq. (fT3|) . 
Fig. [7| shows the "standard Glauber A^coii eccentricity". Another more sophisticated 
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Fig. 6. A comparis on o f the elliptic flow of pions and protons to the elliptic flow of the multi-strange 
(h and f2~ hadronJSI. 



model is provided by the KLN model whi ch is b ased on the ideas of gluo n satur ation 
and the Color Glass Condensate (CGC) 144145] implemented in Refs. SSHZl. This 
model is a safe upper bound on what can be expected for the eccentricity from 
saturation physics and is also shown in Fig. [T] We can not describe the details 
of this model and its implementation here. However, the physical reason why this 
model has a sharper eccentricity is the readily understood: the center of one nucleus 
(nucleus A) passes through the edge of the other nucleus (nucleus B). Since the 
density of gluons per unit area in the initial wave function is larger in the center 
of a nucleus relative to the edge, the typical momentum scale of nucleus A (^ 
Qs,a) is larger nucleus B (~ Qs,b)- It is then difhcult for the long wavelength (low 
momentum) gluons in B to liberate the short wavelength gluons in A. The result is 
that the production of gluons falls off more quickly near the x edge relative to the y 
edge making the eccentricity larger. Clearly this physics is largely correct although 
the magnitude of the effect is uncertain. Another CGC estimate of e is based on 
classical simulations of Yang-Mills fields. The simulations include similar saturation 
physics but model the production and non-perturbative sectors differently. The 
eccentricity from these simulations is also shown in Fig. [7] and is similar to the 
Nco\i eccentricity °. Thus the predictions of the KLN model seem to be a safe 
upper bound for the eccentricity in heavy ion collisions. Note that an important 
phenomenological consequence of the the KLN model is that the eccentricity grows 
with beam energy and is expected to increase about 20% from the RHIC to the 
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b [fm] 

Fig. 7. Figure from Ref.llS! showing various estimates for the initial eccentricity in heavy ion 
coUisions. The physics of the KLN eccentricity is described in the text. In the KLNmodel the 
eccentricity is expected to increase by about 20% when going from RHIC to the LHcJiD. 




Another important aspect in heavy ion coUisions when interpreting the eUip- 
tic flow data is fluctuations in the initial eccentricity. These fluctuations a re not 
accounted for in Fig. [71 The history is complicated and is reviewed in Refs.S^^^. 
There are fluctuations in the initial eccentricity of the participants especially in 
peripheral AuAu and CuCu collisions. Thus rather than using the continuum ap- 
proximation given in Eq. (|13p it is better to implement a Monte-Carlo Glauber 
calculation and estimate the eccentricity using the "participant plane eccentricity" . 
Fig. [8] illustrates the issue: In a given event the ellipse is tilted and the eccentricity 
depends on the distribution of participants. This event by event eccentricity is de- 
noted tpp in the literature. Clearly the experimental goal is to extract the response 
coefficient C relating the elliptic flow to the eccentricity on an event by event basis 

V2 - Ctpp . (18) 

If the flow methods measured (W2), then we could simply divide the measured 
flow to determine the response coefficient, C = (^2) / {tpp). The PHOBOS collab- 
oration deciphered the confusing CuCu data by recognizing the need for epp and 
following this procedurJ^. However, it was generally realized (see in particular. 
Ref.'-' ) that the elliptic flow methods do not measure precisely (W2). Some meth- 
ods (such as two particle correlations V2 {2}) are sensitive to V^wf), while other 
methods (such as the event plane method V2 {EP}) measure something closer to 
(W2). What precisely the event plane method measures depends on the reaction 
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plane resolution in a known way^. So just dividing the measured flow by the av- 
erage participant eccentricity is n ot entirely correct. The appropriate quantity to 
divide by depends on the method I 52 | 50 | 53 [ Gaussian approximation for the 
eccentricity fluctuations this can be worked out analytically. For instance, the two 
particle correlation eUiptic flow V2 {2} (which measures a/ (wf)), should be divided 
by \/ {^"pp) ■ An important corollary of this analysis is that V2 {4} {v2 measured 
from four particle correlations) can be divided by eg of Eq. to yield a good 
estimate of the coefficient C. This is the policy adopted in Fig. [5l Unfortunately, in 
the most peripheral AuAu bins and in CuCu the Gau ssian approximation is poor 
due to strong correlations amongst the participantJ^SI. These correlations arise be- 
cause participants come in pairs and every participant is associated with another 
participant in the other nucleus. Presumably the last centrality bin in Fig. [5] could 
be moved up or down somewhat due to non-Gaussian correct ions of this sort. With 
a complete understanding of what each method measures, Ref.^^was able to make a 
simple model for the fluctuations and non-flow and show that {V2) measured by the 
different methods are compatible to an extremely good precision. This work should 
be extended to the CuCu system where non-G aussian fluctuations are stronger 
and ultimately corroborate the PHOBOS 

analysisPMl. 

This is a worthwhile goal 
because it will clarify the transition into the hydrodynamic regime^. 



2.4. Summary 

In this section we have gone into considerable experimental detail - perhaps more 
than necessary to explain the basic ideas. The reason for this lengthy summary is 
because the trends seen in Fig. [5] were not always so transparent. The relatively co- 
herent hydrodynamic and kinetic interpretation of the observed elliptic flow (which 
was previewed in Section 12.21 and which is discussed more completely below) is the 
result of careful experimental analysis. 
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Fig. 9. Figure from Ref.ESl illustrating the energy density and pressure by of QCD computed 
with A^r = 8 lattice data. (In this figure e is energy density e(T) and the pressure p is denoted 
with V throughout this review.) egs/T^ = eg^/T* is the energy density of a free three flavor 
massless QGP (see text). 

3. The Shear Viscosity in QCD 

In this section we will discuss thermal QCD in equilibrium with the primary goal 
of collecting various theoretical estimates for the shear viscosity in QCD. 

The prominent feature of QCD at finite temperature is the presence of an ap- 
proximate phase transition from hadrons to quarks and gluons. The Equation of 
State (EoS) from lattice QCD calculations is shown in Fig. [9l and the energy den- 
sity e(T) shows a rapid change for the temperature range, T ~ 170 — 220 MeV. As 
estimated in Section HI the transition region is directly probed during high energy 
heavy ion collisions. 

Well below the phase transition, the gas of hadrons is very dilute and the ther- 
modynamics is dominated by the measured particle spectrum. For instance the 
number of pions in this low temperatures regime is 

/■ d3p 1 

~ " 7 (2^)3 e^^/T - 1 ' ^ ' 

where Ep ~ ^Jp^ -t- and (i,r = 3 counts the three fold isospin degeneracy, 
TT^, TT", tt'^, in the spectrum. If all known particles are included up to a mass 
mros < 2.5 GeV, the resulting Hadron Resonance Gas (HRG) equation of state does 
a reasonable job of reproducing the thermodynamics up to about T ~ 180 MeV. 
Howe ver, the validity of this quasi-particle description is unclear above a tempera- 
ture of^, T ~ 140 MeV. As the temperature increases, the hadron wave functions 
overlap until the medium reorganizes into quark and gluon degrees of freedom. Well 
above the transition the QCD medium evolves to a phase of massless quarks and 
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gluons. The energy density is approximately described by the Stefan-Boltzmann 
equation of state 

- flgluc J ^2^)3 gfip/T _ 1 ' ^q""k - flquark J ^^^y ^^^/j, _^ ^ , ) 

where dgiuo = 2x8 counts spin and color, and dquark = 2x2x3x3 counts spin, anti- 
quarks, flavor, and color. Performing these integrals we find, 655 = egiuo + Cquark — 
15.6 as illustrated by the line in the top-right corner of the figure. 

We have described the particle content well above and well below the transition. 
Near the approximate phase transition the validity of such a simple quasi-particle 
description is not clear. The transition is a rapid cross-over where hadron degrees 
of freedom evolve into quark and gluon degrees of freedom rather than a true phase 
transition. All correlators change smoothly, but rapidly, in a temperature range of 
T ~ 170 — 210 MeV. From a phenonienological perspective the smoothness of the 
transition suggests that the change from quarks to hadrons should be thought of 
as a soft process rather than an abrupt change. 

Lattice QCD simulations have determined the equation of state rather well. 
However, in addition to the equation of state, we need to estimate the transport 
coefficients to assess whether the heavy ion reactions produce enough material, over 
a large enough space-time volume to be described in thermodynamic terms. The 
shear and bulk viscosities govern the transport of energy and momentum and are 
clearly the most important. 

Later in Section [4] and Section [5] we will describe the role of shear viscosity in 
the reaction dynamics. In this section we summarize the shear viscosities found 
in various theoretical computations which will place these dynamical conclusions 
in context. A good way to implement this theoretical summary is to form shear 
viscosity to entropy ratid^, 77/s. To motivate this ratio we remark that it seems 
difficult to transport energy faster than a quantum time scale set by the inverse 
temperatur^, 

h 

Equant ^ 7^ • 

kbT 

A sound wave propagating with speed Cs will diffuse (or spread out) due to the shear 
viscosity. Linearized hydrodynamics shows that this process is controlled by the 
momentum diffusion coefficient, = rj/le+V), where e+V is the enthalpy (see for 
example Ref.t^^. Noting that the diffusion coefficient has units of (distance)^ /time, 
a kinetic theory estimate for the diffusion process yields 

„2 



A, = ~ , (21) 



where tu is the particle relaxation time and v'^^i ~ is the particle velocity. Dividing 



by v^Yi and using the thermodynamic estimates 

sT r-. ev^i^^V ^nkeT, (22) 



^In this paragraph we will restore h and the Boltzmann constant, fcs. 



May 14, 2009 21:38 WSPC/INSTRUCTION FILE hydro 



Viscous Hydrodynamics and the Quark Gluon Plasma 15 

we see that 

- ^trT^ ^ . (23) 

S Kb Tquant 

Therefore, 77/5 is the ratio between the medium relaxation time and the quantum 
time scale Tquant in units of h/ks, i.e. a measure of the transport time in "natural 
units" . 

In the dilute regime the ratio between the medium relaxation time and the 
quantum time scale is long and kinetic theory can be used to calculate the shear 
viscosity to entropy ratio. First we consider a simple classical massless gas with 
particle density n and a constant hard sphere cross section ctq. The equation of 
state of this gas is e = SP ~ ZnT and the shear viscosity is computed using kinetic 
theorjl^ 

7,-1.2-, (24) 

<Jo 

The entropy is s = (e + "P) /T and the resulting shear to entropy ratio is 

^-0.3 — . (25) 
s noo 

In what follows, this calculation will provide a qualitative understanding of more 
sophisticated kinetic calculations. 

In the dilute hadronic regime, 77/5 was calculated in using measured 

elastic cross sections for a gas of pions and kaons. In the tttt phase shifts there is a 
prominent p resonance, while in the t:K channel there is a prominent K* resonance. 
Thus the equation of state of this gas is well modeled by an ideal gas of tt, p and 
K* '^^^ The viscosity of this mixture was computed in Ref.''^' and the current 
author digitized this viscosity, computed the entropy, and determined the rjj s ratio. 
This is shown in Fig. [TOl Slightly larger values were obtained in Rcf."^^ which also 
estimated the range of validity for hadronic kinetic theory, T < 140 MeV. Finally 
a more involved Kubo analysis of the UrQMD hadronic transport model 1^ (which 
includes many resonances) is also displayed in Fig. 1101 

At asymptotically high temperatures the coupling constant is weak and the 
shear viscosity can be computed using perturbation theory. Initially, only 2^2 
elastic scattering was considered, and the shear viscosity was computed ina lead- 
ing log plasma with self consistent screeninj^. Later it was recognizecP^^^ that 
coUinear Bremsstrahlung processes are important for the calculation of shear viscos- 
ity and this realization ultimately resulted in a complete leading order calculation 

We can estimate rj/s in the perturbative plasma using Eq. ([24| with s (x 
and a oc aj./T'^, 

(26) 

The final result from a complete calculation has the form 

!l^J_Fimn/T), (27) 
s at 
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Fig . 10 . (Color Online) A compilation of values of rj/s. The results from Prakash et al are from 
ReflSSl and describe a meson gas of pions and kaons (and indirectly K* and p) computed with 
measured cross sections. The black poi nts are based on a Kubo analysis of the UrQMD code 
which includes many higher resonances EZl, The red lines are different implementations of the 
AMY (Arnold, Moore, Yaffe) calculation of shear viscosity 1221, In each curve the Debye scale is 
fixed mo = 2T. In the dashed red curves the (one loop three flavor) running coupling is taken at 
the scale /x. In the solid red curves as is kept fixed. The two loop running coupling is shown with 
= 2ttT for comparison and the two loop = ttT (not shown) is similar to the one loop = 2ttT 
result. In the AMY curves, changing the Debye mass by it0.5T changes ri/s by ^ ±30%. Finally 
the thin dashed line indicates a simple model discussed in the text with ^^fp = 1/T. 



where F{mD/T) is a function of the Debye mass which was computed for mu/T 
small and then extrapolated to more realistic values There are many scales in 
the problem and it is difficult to know what precisely to take for the Debye mass 
and the coupling constant. At lowest order in the coupling, the Debye mass iJ^ 



f + ^')5^T^^ (28) 



but this is too large to be considered reliable. For definiteness we have evaluated 
the leading coupling constant in Eq. l[27|) at a scale of ttT and set the Debye mass to 
rriY) — 2T. The resulting value oirj/s is shown in Fig. [TOl Various other alternatives 
are explored in the figure and underscore the ambiguity in these numbers. 

Clearly all of the calculations presented have a great deal of uncertainty around 
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the phase transition region. On the hadronic side there are a large number of inelas- 
tic reactions which become important. On the quark gluon plasma side, the strong 
dependence on the Debye scale and the coupling constant is disconcerting. It is 
very useful to have a strongly coupled theory where the shear viscosity to entropy 
ratio can be computed exactly. In strongly coupled Af = 4 SYM theory with a large 
num ber of c olors , rj/s can be computed using gauge gravity duality and yields the 

resuijmsi 

— f- (29) 

From the perspective of heavy ion physics this result was important because it 
showed that there exist field theories where rj/s can be this low. Although J\f — 
4 has no particle interpretation, we note that extrapolating Eq. (|24p by setting 
^mfp — i/nao = 1/nT yields a value for rj/s which is approximately equal to the 
SYM result. In Fig. [10] we have displayed this numerology with fmtp = 1/T for 
clarity. 

There are many aspects of transport coefficients which have not been reviewed 
here. For instance, ther e is an ongoing effort to determine the transport coefficients 
of QCD from the lattice'^^'^. While a precise determination of the transport coef- 
ficients is very difficult'" '-^ the lattice may be able to determine enough about 
the spectral densities to distinguish the orthogonal pictures represented by A/" = 4 
SYM theory and kinetic theory"". This is clearly an important goal and we refer 
to Ref.l^ for theoretical background. Also throughout this review we have empha- 
sized the shear viscosity and neglected the bulk viscosity. This is because on the 
hadronic side of the phase transition the bulk viscosity is a thousand tim es s maller 
than the shear viscosity in the regime where it can be reliably calculated 1^21. Simi- 
larly on the high temperature QGP side of the phase transition the bulk viscosity 
is also a thousand times smaller than sheaJ^. Ho wever, nea r a second order phase 
transition the bulk viscosity can become very largJ^HZlIZSl, Nevertheless the rapid 
cross-over seen in Fig. [9] is not particularly close to a second order phase transition 
and universality arguments can be questioned (see for a discussion in the 

context of the chiral susceptibility.) Given the ambiguity at this moment it seems 
prudent to leave the bulk viscosity to future review. 



4. Hydrodynamic Description of Heavy Ion Collisions 

In the previous sections we analyzed the phase diagram of QCD and estimated the 
transport coefficients in different phases. In this section we will study the hydrody- 
namic modeling of heavy ion collisions. 

In Section [4.21 we will consider ideal hydrodynamics and assume that the mean 
free paths are small enough to support this interpretation. Subsequently we will 
study viscous hydrodynamics in Section [4?3l Section [4^ will analyze the ratio of the 
viscous terms to the ideal terms and use the estimates of the transport coefficients 
given above to assess the validity of the hydrodynamic interpretation. Section [4.61 
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will discuss the recent advances in interpreting the hydrodynamic equations beyond 
the Navier Stokes limit. This work will lay the foundation for the more detailed 
hydrodynamic models presented in Section [6] 

4.1. Ideal Hydrodynamics 

The stress tensor of an ideal fluid and its equation of motion are simply 

T'"' ^ eu^u" + VA"" , df.T"" = , (30) 

where e is the energy density, 'P(e) is the pressure, and = (7,7V) is the four 
velocity. Here we will use the metric (—,+,+,+) and define the projection tensor, 
A^"^ = 5^" + u^u", with w^M^ = -1 and A^^u^ ^ 0. This decomposition of 
the stress tensor is simply a reflection of the fact that in the local rest frame of 
a thermalized medium the stress tensor must have the form, diag(e, "P, 7^, P). In 
developing viscous hydrodynamics we will define two derivatives which are the time 
derivative D, and the spatial derivatives in the local rest frame 

D EE u^a^ , = A^'-a^ . (31) 

Using — —Ufj_D + and Ufj^Du'^ — 0, the ideal equations of motion can be 
written 

De = -ie + V)V^u^', (32) 

Du^^-—-. (33) 

e + V 

The first equation says that the change in energy density is due to the PdV work 
or equivalently that entropy is conserved. To see this we associate V ^u^ with the 
fractional change in volume per unit time in the co-moving frame, dV/V ^ dt x 
Vfj^u^, and use the thermodynamic identity, d{eV) — Td{sV) — VdV . The second 
equation says that the acceleration is due to the gradients of pressure. The enthalpy 
plays the role of the mass density in a relativistic theory. 



4.2. Ideal Bjorken Evolutions and Three Dimensional Estimates 

In this section we will follow an analysis due to BjorkerP^ and apply ideal hydro- 
dynamics to he avy ion c ollisions. Bjorken's analysis was subsequently extended in 
important wayJ^^IISIlSl in a high energy heavy ion collisions the two nuclei pass 
through each other and the partons are scarcely stopped. This statement under- 
lies much of the interpretation of high energy events and an enormous amount of 
data is consistent with this assumption. For a time which is short compared to the 
transverse size of the nucleus, the transverse expansion can be ignored. 

Given that the nuclear constituents pass through each other, the longitudinal 
momentum is much much larger than the transverse momentum. Because of this 
scale separation there is a strong identification between the space-time coordinates 
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and the typical z momentum. For example a particle with typical momentum 
and energy E will be found in a definite region of space time 

v^-^?^^'-. (34) 
E t ^ ' 

This kinematics is best analyzed with proper time and space-time rapidity vari- 
able^, T and rjs 

1, ft + z 



2 °\t~ z 

At a proper time r particles with rapidity y are predominantly located at space 
time rapidity 77^ 

1 p,+E 1 t + z 

y^o^og ~ - log EE ris . 35) 

2 E — pz 2 t — z 

Fig. [TT] illustrates these coordinates and shows schematically the identification be- 
tween rjg and y. At an initial proper time Tq, there is a collection of particles 
predominantly moving with four velocity in each space-time rapidity slice 

1 



The beam rapidity at RHIC is j/bcam — 5.3 and therefore roughly speaking the 
particles are produced in the space-time rapidity range —5.3 < rjs < 5.3. It is 
important to realize that (up to about a unit or so) each space-time rapidity slice 
is associated with a definite angle in the detector. For ultra-relativistic particles 
E p we have 

1, fp + pA 1, /l + cos0\ 
rys ~ y ~ - log = o log "j n = '?p=-i= ' (37) 



2 \P — Pz J 2 \1 — cos( 

where a particular 6 is shown in Fig. [1] The measured pseudo-rapidity distribution 
of charged particles is shown in Fig. [121 We can estimate the energy in a unit 
of pseudo-rapidity by taking (E) ~ 0.5 GeV as the energy per particle. Then the 
energy in a pseudo-rapidity unit is 

- (E) /^"^ X 1.5 ~ 3.0 GeV x {Np/2) , 



where {Np/2) ~ 170 is the number of participant pairs in a central event. The 
factor of 1.5 has been inserted to account for the fact that there are approximately 
equal numbers of tt"*" , tt" and 7r° (the most abundant particle) but only 7t~^ and 7r~ 
are counted in diVch/df/paoudo- This estimate agrees reasonably with the measured 
dEr/dri^,,^^, ~ 3.2 GeV x Np/2 from Ref.lSHl. 



^ Here r}s denotes the space time rapidity, Jjpaeudo denotes the pseudo-rapidity (see below), r] 
denotes the shear viscosity. In raised space time indices in t, r]s coordinates we will omit the "s" 
when confusion can not arise, e.g. -k^^ = n^'^' . 
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Fig. 11. A figure motivating for the Bjorken model. The space between the dashed fines of constant 
rjs are referred to as a space-time rapidity slice in the text. Lines of constant proper time t are 
given by the solid hyperbolas. The collection of particles in the rjs = rapidity slice is indicated 
by the small arrows for the central (tJs = 0) rapidity slice only. The solid arrows indicates the 
average four velocity in each slice. The spectators are those nucleons which do not participate 
in the collision and lie along the light cone. 
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Fig. 12. The measured charged particle pseudo-rapidity distribution dA'^h/'^'JpEioudo for different 
beam energies divided by the number of participant pairs, Np/2. Np/2 ~ 170 for a central (0-6%) 
AuAu collision. This review focuses on y/s = 200 GeV/nucleon. 



Bjorken used these kinematic ideas to estimate the initial energy density in the 
r]s = rapidity shoe at an initial time, Tq ~ 1 fm. The estimate is based on the fairly 
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well supported assumption that the energy which finally flows into the detector 
di?T / drj-p^^^^a largely reflects the initial energy in a given space-time rapidity slice 
_ 1 A£' _ 1 A£; _ 1 dEr 

^5.5^. (39) 
tm 

In the last line we have estimated the area of a gold nucleus as A ~ 100 fm^, taken 
To — 1 fm, and used the measured cIEt / drj^^^^^^ . This estimate is generally consid- 
ered a lower limit since during the expansion there is VdV work as t he particl es in 
one rapidity slice push against the particles in another rapidity slicJi^^SlZSl (|See 
Fig. [TT|) . Using the equation of state in Fig. [9] we estimate an initial temperature, 
T{to) — 250 MeV. As mentioned above this estimate is somewhat low for hydro- 
dynamic calculations and a more typ ical temperature is T ~ 310 MeV, which has 
roughly twice the Bjorken density24l 

As seen in Fig. [121 the distribution of the energy density e(To,r?s) in space- 
time rapidity is not uniform. In the Color Glass Condensate (CCC) picture for 
instance, the final distri buti on of multiplicity is related to the x distribution of 
partons inside the nucleuJ^H. Bjorken made the additional simplifying assumption 
that the energy density is uniform in space-time rapidity, i.e. e{To, r]s) — siTo)- With 
this simplification, the identification between the fluid and space time rapidities 
remains fixed as the fiuid flows into the forward light cone. 

We have discussed the motivation for the Bjorken model. Formally the model 
consists of the following ansatz for the hydrodynamic variables 

e{t, x) = e(r) , x) = u^, u^', u^) = (cosh(77,), 0, 0, sinh(77,)) . (40) 

The model is invariant under boosts in the z direction. Thus given a physical 
quantity at mid-rapidity (77^ = 0), one can determine this quantity at all other 
rapidities by a longitudinal boost. We will use curvilinear coordinates wherJSl 

x'^ = (T,Xi,77,), .g^, = diag(-l,l,l,r2), (^i^ u^ u^ u") = (1, 0, 0, 0) . 

(41) 

In this coordinate system boost invariance implies that everything is indepen- 
dent of rjg. To interpret a tensorial component in these coordinates, we multiply 
by ^/9rrn — for every raised ry^ index, and subsequently associate the prod- 
uct with the corresponding cartesian component at mid-rapidity. For example, 
^221)))) _ T^^l^^^p = 7?. Similarly, rw'' u^\r^,=o ~ for boost invariant flow. 

Substituting the boost invariant ansatz (Eq. (|40p ) into the conservation laws 
yields the following equation for the energy densitjo 

^^^-^^ (42) 

dr r 

'^A quick way to derive this is to work in a neighborhood of z = rjg = where ~ z/t. Substituting 
this approximate form into df^Tf^" = in cartesian coordinates, quickly yields Eq. I I42I I with the 
replacement t ^ t. 
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Multiplying this equation by the volume of a space time rapidity slice V = TArjsA 
(see Fig. [11]) we find 

d(e TArjsA) = -V d(rA?7, A) , (43) 

and we can interpret this result - as saying that the energy energy per unit space- 
time rapidity {erArisA) decreases due to the VdV work. It is for this reason that 
the Bjorken estimate esj (which assumes that the r.h.s. of Eq. equals zero) 
should be considered a lower bound. In general, assuming Bjorken scaling (Eq. (j40p ) 
and the conservation laws, but not assuming local thermal equilibrium, one finds 

de e + T^^ , , 

■X- = . (44) 

dr r 

where T^^ = t^T^^ is the effective longitudinal pressure. Viscous corrections will 
modify T^^ from its equilibrium value of V. 

Returning to the equilibrium case, Eq. (|42| can be solved for a massless ideal 
gas equation of state (e = 3V oc T*) and the time dependence of the temperature 
is 

/r \ 1/3 

nr)^n{^) , (45) 

where To is the initial temperature. The temperature decreases rather slowly as a 
function of proper time during the initial one dimensional expansion. This will turn 
out to be important when discussing equilibration. For a massless ideal gas, the 
entropy is s — {e + V) /T oc and decreases as 

sir) ^ So-. (46) 
r 

Now we discuss what happens when the initial energy density distribution is not 
uniform in rapidity. Due to pressure gradients in the longitudinal direction, there 
is some longitudinal acceleration. This changes the strict identification between the 
space time rapidity and the fluid rapidity given in Eq. (l36l) . It also changes the 
temperature dependence given above. One way to quantify this effect is to look at 
the results of 3D ideal hydrodynamic calculations and study the differences between 
the initial energy distribution in space-time rapidity J d^xj_ e(To, x_l, 77s) and the 
final energy distribution, J d^x^ e(r/, x_l, rys)- Generally, the final distribu tion in 
space-time rapidity is similar to the initial distribution in space-time rapidity ^1^^^. 
Therefore, the effect of longitudinal acceleration is unimportant until late times. 

The nuclei have a finite transverse size, Rau 6 fm. After a time of order 

r , 

c 

the expansion becomes three dimensional. To estimate how the temperature evolves 
during the course of the resulting 3D expansion, consider a sphere of radius R which 
expands in all three directions. The radius and volume increase as 

R(XT, V(XT^. 
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Since for an ideal expansion the total entropy in the sphere is constant, the entropy 
density decreases as 1/r'^ and the temperature decreases as 

11 

s oc — , r cx — . (47) 

T 

Here we have estimated how the entropy decreases during a one and three 
dimensional expansion of an ideal massless gas. Now if during the course of the 
collision there are non-equilibrium processes which generate entropy that ultimately 
equilibrates, the temperature of this final equilibrated gas will be larger than if the 
expansion was isentropic. Effectively the temperature will decrease more slowly. To 
estimate this effect in a one dimensional expansion, we imagine a free streaming 
gas where the longitudinal pressure is zero. Then from Eq. ()42|) we have 

(48) 

In the sense discussed above, this equation may be integrated to estimate that the 
temperature and entropy decrease as 

T cx -ij , s cx ^ . (49) 

Similarly in a three dimensional expansion we can estimate how entropy production 
will change the powers given in Eq. (jiTj) . Again consider a sphere of radius R which 
expands in all three directions, such that R oc t and ^ cx r^. For a free expansion 
without pressure the total energy in the sphere is constant, and the energy density 
decreases as Similarly, we estimate that the temperature and entropy density 

decrease as 

T cx , 5 oc . (50) 



In summary we have estimated how the temperature and entropy density de- 
pend on the proper time r during the course of an ideal and non-ideal ID and 31? 
expansion. This information is recorded in Table [T] These estimates are also nicely 
realized in actual hydrodynamic simulations. Fig. [13] shows the dependence of en- 
tropy density as a function of proper time r. The figure indicates that the entropy 
decreases as 1 /r during an initial one dimensional expansion and subsequently de- 
creases as when the expansion becomes three dimensional at a time of order 
^ 5 fm. These basic rules will be useful when estimating the relative size of viscous 
terms in what follows. 



4.3. Viscous Bjorken Evolution and Three Dimensional Estimates 

This section will analyze viscosity in the context of the Bjorken model with the 
primary goal of assessing the validity of hydrodynamics in heavy ion collisions. In 
viscous hydrodynamics the stress tensor is expanded in all possible gradients. Using 
lower order equations of motion any time derivatives of conserved quantities can 
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Fig. 13. Figure from Kef.! showing the entropy density (s) in CuCu simulations as a function 
of proper time r using ideal and viscous hydrodynamics. The top set of lines shows the entropy 
in the center of the nucleus-nucleus collision, (r = Ofm), and the bottom set of lines shows the 
analogous curves closer to the edge (r = 3fm). During an initial one dimensional expansion the 
entropy density decreases as s oc l/r. Subsequently the entropy decreases as s oc l/r^ when the 
expansion becomes three dimensional at a time, r ^ 5 fm. The lines labeled by (0 + 1) ideal and 
(0 + 1) viscous are representative of the ideal and viscous Bjorkcn results Eq. I I42I I and Eq. II55I I 
respectively. 



Quantity 


ID Expansion 


3D Expansion 


T 






s (X 







Table 1. Dependence of temperature and entropy as a function of time in a ID and 3D expansion. 
The indicated range, for instance 1/3 1/4, is an estimate of how extreme non-equilibrium effects 
could modify the ideal power from 1/3 to 1/4. 



be rewritten as spatial derivatives. First the stress tensor is decomposed into ideal 
and viscous pieces 

T^" = T;^",! + TT^" + UAf"" , (51) 

where Tj^'' is the ideal stress tensor (Eq. ((30|) ) and 11 is the bulk stress, tt^'' is 
the symmetric traceless shear tensor and satisfies the orthogonality constraint, 
■K^'^Uy = 0. The equations of motion are the conservation laws dfj.T^" = to- 
gether with a constituent relation. The constituent relation expands tt^'' and 11 
in terms gradients of the conserved charges T^" and T'^* or their thermodynamic 
conjugates, temperature T and four velocity . To first order in this expansion. 
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the equations of motion are 

a^T^'^^o, TT^"^ = -7^(7^% n = -cv^u^, (52) 

where 77 and C, are the shear and bulk viscosities respectively, and we have defined 
the symmetric traceless combination 

cr^"^ = V^'ii'^ + - \l\^"'Vxu^ . (53) 

For later use we also define the bracket (. . .) operation 

{An = Ia'^-A''^ {A^p + A;,^) - Ia^^A^'^A^;, , (54) 

which takes a tensor and renders it symmetric, traceless and orthogonal to . Note 
that cr'"' = 2{d^'v}'). 

We now extend the Bjorken model to the viscous case following ReffEl. The 
bulk viscosity is neglected in the following analysis and we refer to Section [3] for 
a more complete discussion. Substituting the Bjorken ansatz (Eq. (j40|) ) into the 
conservation laws and the associated constituent relation (Eq. (j52|) ) yields the time 
evolution of the energy density 

de^_e^rV-y^_ (55) 
ctr T 

The system is expanding in the z direction and consequently the pressure in the z 
direction is reduced from its ideal value. Formally this arises due to the gradient 
dzU^ = 1/t and the constituent relation Eq. ([52]) 

T^^ = 7?_l!^. (56) 

3 T 

Thus during a viscous Bjorken expansion the system will do less longitudinal work 
than in the ideal case. 



4.4. The Applicability of Hydrodynamics and ri/s 

Comparing the viscous equation of motion Eq. ([55)1 to the ideal equation of motion 
Eq. ([m, we see that the hydrodynamic expansion is controlled by 

:-«l. (57) 



e + Vr 

This is a very general result and is a function of time and temperature. Using the 
thermodynamic relation e + V ^ sT, we divide this condition into a constraint on 
a medium parameter rj/s and a constraint on an experimental parameter 1/tT 

- X «1. (58) 

medium parameter experimental parameter 

If the experimental conditions are favorable enough, it is appropriate to apply 
hydrodynamics regardless of the value of 77/s. This is the case for sound waves in air 
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where although rj/s is significantly larger than the quantum bound, hydrodynamics 
remains a good effective theory. However, for the application to heavy ion collisions, 
the experimental conditions are so unfavorable that only if rj/s is close to the 
quantum bound will hydrodynamics be an appropriate description. 

For instance, we estimated the experimental condition in Section 14.11 

1 /lfm\ /300MeV\ , , 

Here we have evaluated this experimental parameter at a specific initial time Tq and 
will return to the time evolution of these estimates in the next section. In Section [3] 
we estimated the medium parameter rj/s and can now place these results in context 

From this condition we see that hydrodynamics will begin to be a good approxima- 
tion for r]/s ^ 0.3 or so. This estimate is borne out by the more detailed calculations 
presented in Section [S] Reexamining Fig.fTUl we see that the value of 77/ s ~ 0.3 is at 
the low end of the perturbative QGP estimates given in the figure and it is difficult 
to reconcile the observation of strong collective flow with a quasi-particle picture of 
quarks and gluons. Thus the estimates of rj/s coming from the RHIC experiments, 
which are based on the hydrodynamic interpretation of the observed flow, should 
be accepted only with considerable care. 



4.5. Time Evolution 

In the previous section we have estimated the applicability of hydrodynamics at a 
time To « 1 fm. In this section we will estimate how the size of the viscous terms 
depends on time. For this purpose we will keep in mind a kinetic theory estimate 
for the shear viscosity 



and estimate how the gradient expansion parameter in Eq. (|55|1 depends on time. 
We will contrast a conformal gas with a cx 1/T^ (e-5- perturbation theory or A/" = 4 
SYM) to a gas with fixed cross section, a = ao- There are clearly important scales 
in the quark gluon plasma as the medium approaches the transition point. For 
instance spectral densities of current-current correlators near the transition point 
show a very d iscernible correlation where the p meson will form in the hadron 
p]^g^gg[83l84l69|^ Thus the intent of studying this extreme limit with a constant cross 
section is to show some of the possible effects of these scales. Further the constant 
cross section kinetic theory has been used to analyze the centrality dependence of 
elliptic 

First consider a theory where the temperature T is the only scale and also 
consider a ID Bjorken expansion. The shear viscosity is proportional to and the 
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enthalpy scales as T**, so the hydrodynamic expansion parameter scales as 

(62) 

{e + V)T tT t2/3 ■ ^ 

In the last step we have used the fact that for a scale invariant gas undergoing 
an ideal Bjorken expansion the temperature decreases as l/r^^'^. In general if we 
have some non-equilibrium processes which produce entropy during the course of 
the expansion, the temperature will decrease more slowly than estimated in the 
ideal gas case - see Table [TJ The result is that we do not expect the temperature 
to decrease more slowly than 1/t^/'*, and we can estimate that the hydrodynamic 
expansion parameter evolves as 

^ « « • (63) 



(e + -P) r tT t2/3-3/4 

Thus during a ID expansion of a conformal gas the system will move closer to 
equilibrium. 

Compare this scale invariant theory to a gas with a very definite cross section 
Uo- For a constant Uq the hydrodynamic expansion parameter evolves as 

1 1 /I \ 0+1/4 

^ loc^ocfi) . (64) 



[e + V) T sar 

Thus, with a constant cross section, the gas will move neither away nor toward 
equilibrium as a function of time. Non-equilibrium physics will make the matter 
evolve slowly toward equilibrium. 

Now we will compute the analogous effects for a three dimensional expansion. 
In the conformal case r] oc and T oc ^ , so that the final result is 

1 1 /l\"^^^'' 
oc — cx - ) . (65) 



{e + V)T tT 

Thus a conformal gas expanding isentropically in three dimensions also moves nei- 
ther away nor towards equilibrium, though entropy production will cause it to slowly 
equilibrate. Similarly for gas with a constant cross section the hydrodynamic pa- 
rameter evolves as 



{e + V)T S(JT ^ ' 

In estimating this last line we have used Table [TJ Thus we see that a gas with fixed 
cross-sections which expands in three dimensions very rapidly breaks up. 

The preceding results are summarized in Table [51 Essentially the heavy ion colli- 
sion proceeds along the following line of reasoning. First, there is a one dimensional 
expansion where the temperature is the dominant scale in the problem. The param- 
eter which controls the applicability of hydrodynamics 'rj/[{e -\- 'P)t\ decreases as a 
function of time; hydrodynamics gets better and better, evolving according to the 
upper left corner of Table [2] As the system expands and cools toward the transition 
region additional scales enter the problem. Typically at this point r ^ 4fm/c the 



May 14, 2009 21:38 WSPC/INSTRUCTION FILE hydro 



28 Derek A. Teaney 



Model 



ID Expansion 



3D Expansion 



rj (X 





Table 2. Dependence of the hydrodynamic expansion parameter r;/[(e + 'P)t] as a function of 
time for two different functional forms for rj {rj oc T and rj oc T^) and two expansion types (ID 
and 3D). A range of powers is given; the first power corresponds to ideal hydrodynamics and the 
second power corresponds to an estimate of non-equilibrium evolution. 

expansion also becomes three dimensional. The system then enters the lower right 
corner of Table [2] and very quickly the nucleus- nucleus collision starts to break up. 
We note that it is necessary to introduce some scale into the problem in order to see 
this freezeout process. For a conformal liquid with ij (x the system never freezes 
out even for a 3D expansion. This can be seen by looking at the upper-right corner 
of the table and noting that the hydrodynamic expansion parameter behaves as 



and therefore approaches a constant (or slowly equilibrates) at late times. From 
this discussion we see that the temperature dependence of the shear viscosity is 
ultimately responsible for setting the duration of the hydrodynamic expansion. 

4.6. Second Order Hydrodynamics 

In the previous sections we developed the first order theory of relativistic viscous 
hydrodynamics. In the first-order theory there are reported instabilities which are 
associated with the gradient expansion^. Specifically, in the first-order theory the 
stress tensor is instantly specified by the constituent relation and this leads to 
acausal propagatiorP^ and ultimately the instability. Nevertheless, it was generally 
understood that one could write down any relaxation model which conserved en- 
ergy and momentum and which included some notion of entropy, and the results of 
such a model would be indistinguishable from the Navicr Stokes cquationJSIIHHISS]^ 
Many hydrodynamic models were written down^ ^ ^'^ ^"^ s tarting with a phe- 
nomenological model by Israel and StewarlPSEOI ^nd MiilleJ^. For example in the 
authors own work the strategy was to write down a fiuid model (based on Ref.^^Sf) 
which relaxed on some time scale to the Navier Stokes equations, solve these model 
equations on the computer, and finally to verify that the results are independent of 
the details of the model- Thus the goal was to solve the Navier Stokes equations 
and to estimate the effects of higher order terms. 

Recently an important work by R. Baier, P. Romatschke, D. T. Son, 
A. O. Starinets and M. A. Stephanov (hereafter BRSSS) clarified and classified 
the nature of these higher order terr nj^^J An important impetus for this work came 




(67) 



from the AdS/CFT corres; 




'. Many of the fluid models discussed 
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above were motivated by kinetic theory. However, in the strongly coupled M — A 
plasma, kinetic theory is not applicable, and the precise meaning of these models 
was vague. BRSSS determined precisely in what sense these second order viscous 
equations are theories and in what sense they are models. Simultaneously the Tatta 
group completed the calculation of the second order transport coefficients in A/" = 4 
SYM theory and clarified the hydrodynamic nature of black branes in the process^. 
The spirit of the BRSSS analysis is the following: 

(1) Write the stress tensor as an expansion in all possible second order gradients 
of conserved charges and external fields which are allowed by the symmetries. 
The transport coefficients are the coefficients of this gradient expansion. 

(2) In this expansion temporal derivatives can be rewritten as spatial derivatives 
using lower order equations of motion. 

(3) The conservation laws d^T^'^ = and the associated constituent relation dic- 
tates the dynamics of the conserved charges in the presence of the external field. 
By adjusting the transport coefRcients, this dynamics will be able to reproduce 
all the retarded correlators of the microscopic theory. 

In general, for a theory with conserved baryon number there are many terms. By 
focusing on a theory without baryon number and also assuming that the fluid is 
conformally invariant, the number of possible second order terms is relatively small. 
The classification of gradients in terms of their conformal transformation properties 
was very useful, both theoretically and phenomenologically. At a theoretical level 
there are a manageable number of terms to write down. At a phenomenological 
level the gradient expansion converges more rapidly when only those second or- 
der terms which are allowed by conformal invariance are included (see Section 16]) . 
Subsequently when additional conformal breaking terms are added, the conformal 
classification provides a useful estimate for the size of these terms, i.e. quantities 
that scale as = e — 3p should be estimated differently than those that scale as 
energy density itself. In retrospect, this classification is an "obvious" generalization 
of the first order Navier-Stokes equations. 

Proceeding more technically, in analogy to the constituent relation of the Navier- 
Stokes theory Eq. ([52|) . BRSSS determine that the possible forms of the gradient 
expansion in a conformal liquid are 



+\Aa\a''^) + \2U\n"^) + \Jn\n-^) , (68) 



where the vorticity tensor is defined as 



- ^A'^^A'-'^ {d^up - dpuo.) , (69) 



and c? = 4 is the number of space-time dimensions. Conformal invariance forces a 
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particular combination of second derivatives to have a single coefficient 

(70) 



a — 1 



The time derivative Da'^'^ may be expanded out using lower order equations of 
motion if desired. The constituent relation (Eq. (|68p ) and the conservation laws 
form the second order equations of motion of a conformal fluid. They are precisely 
analogous to the first order theory. As in the first order case, these equations are 
also acausal. 

To circumvent this issue, BRS SS (following the spirit of earlier work by Israel 
and Stewart and Miiller ISH) promote the constituent relation to a dynamical 

equation for the viscous components of the stress tensor tt^'' . Using the lower order 
relation tt^'^ = — rycr^'^, the (conformal) dependence of rj on temperature rj oc T'^~^, 
and the ideal equation of motion Eq. ((32l) . the following equation arises for tt'^'^ 

TT — —rja — 



a — 1 



) + X3 . (71) 



From a numerical perspective the resulting equation of motion is now first order 
in time derivatives, hyper bolic and causal. The modes in this (and similar) models 
have been studied in Refsl ^^ ' ^^ ' ^^ l 

Nevertheless it should be emphasized that the domain of validity of the resulting 
equations is still the same as Eq. (|68p. i.e. the hydrodynamic regime. Thus for 
instance the second order equations should be used in a regime where 

Outside of this regime there is no guarantee that entropy produc tion predicted 
by this model will be positive during the course of the evolution It should 
also be emphasized that this is not a unique way to construct a hydrodynamic 
model which reduces to Eq. ((68|) in the long wavelength limit - see Ref.'^ for an 
example discussed in these terms. What is gua rantee d is that any conformal model 
or dynami cs (suc h as conformal kinetic theory l ^^ l ^^ l or the dynamics predicted by 
AdS/CFT I^^ESl) T^iW be expressible in the long wavelength limit in terms of the 
gradient expansion given above. 

Ther e is an important distinction between the first and second order theories 

| 87 | 72 | 88 | jj^ ^Yie first order theory, the ideal motion is damped, and there are cor- 
rections to the ideal motion of order the inverse Reynolds number 

where L is the characteristic spatial dimension of the system. At is the time of 
observation, and uth is a typical quasi-particle velocity. Thus for sufficiently long 
times the viscous corrections become large and must be resummed by solving the 
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Navier-Stokes equations to capture the damping of the fluid motion. Once this is 
done however, the remaining higher order terms (which are captured by the second 
order theory) are uniformly smafl and modify the Navier-Stokes solution by an 
amount of order 

Often this makes these higher order terms difficult to measure in normal laboratory 
liquids 

For completeness we record th e model equations which have been discussed in 
the heavy ion literature I 5|30 | 9 | 29 | 

(1) The first of these is the simplified Israel-Stuart equation, 

TT^"" = -W'"" - T-n {DtT^'") . (73) 

Since the derivatives do not appear as the combination 

(Dtt^'') + — ^TT^-'V • u , (74) 
a — I 

but rather involve {Dir'^'^) separately, this model does not respect conformal 
invariance. 

(2) The second model is the full Israel-Stewart equation which has the following 
fornP 

= -r^a^'- - {Dnn + \i^^'''^dp (^"'') + 2r. i:''^''^f\ , (75) 



a — 1 



2t,^"(^17'^)^. (76) 



In the last line we have used the conformal relation, tiT/tt^ cx T'^^^ and equation 
of motion, _D(lnT) = —l/{d— 1)V • u. The model is equivalent to taking 
Ai = A3 = and A2 = —'irjTT^ . 

There has been some effort to compute the coefficients of the gradient expansion 
both at strong and weak coupling. The gradient expansion in Eq. (|68p implies that 
the relative size of the coefficients is {imfp/L)'^, and th is is of order [77/(6 -t- V)]^ 
in a relativistic theory. The strong coupling resultj™ are listed in Table 14.61 At 
weak coupling, the results of kinetic calculations are also listed in the table^^ (see 



also Ref .1^3) , Jn kinetic theory the physics of these higher order terms stems from 
the streaming terms and the collision integralJ^SESl To first order in the gradients, 
the distribution is modified from its equilibrium form 

n^n + 6f , (77) 

where 6f oc p^p^aij - see Section [5l Substituting this correction back into the 
Boltzmann equation 



P^d^f = -C[f] , 



(78) 
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Quality 


TV = 4 SYM 


QCD Kinetic Theory 


Relaxation Time 




4-2 ln(2) ~ 2.61 


5.9 to 5.0 (due to g) 


6 


Ai 


2 


5.2 to 4.1 (due to g) 


6 (= V^Tv) 


A2 


-4 ln(2) ~ -2.77 


-11.8 to -10 (= -2?7r^) 


-12 (= -2r]T^) 


As 












Table 3. Compilation of values of (rescaled) second order transport quantities (»?T,r, Ai, A2, A3). All 
numbers in this table should be mul tiplied b y 7?^ / (e -|- V). The complete strong coupling res ults 
are from an amalgamation of Ref.1251 and RefESI, The weak coupling results are from Ref .[2S1 and 
the relaxation time approximation was studied in 

Ref.|25| 

and clarified in Ref.|26|. Hydrodynamic 
simulations of the heavy ion event are not sensitive to these values. In a theory where Ai = rjT-^ 
the second order corrections to a viscous 4-1 dimensional Bjorken evolution vanish. 



leads to several terms which are responsible for the second order corrections to 
hydrodynamics. We enumerate these contributions: 

(1) The T,r and A2 terms are the result of streaming of the first viscous correction 
P^dfiSf and do not involve the collision integral. The common origin of these 
terms ultimately explains the relation between them, A2 — —2'qTT;. 

(2) The contribution to Ai (the visco-elastic tttt term) reflects the streaming 
P^d^Sf and the non-linearities of the collision integral, C[5f]. 

(3) Finally the vorticity-vorticity term does not appear on the LHS of the Bq ltz- 
mann equation and therefore this term vanishes in kinetic theory | 25 | 27|26 | jj^ 
the strong coupling limit the absence of a vorticity-vorticity coupling is not 
understood. 

In the relaxation time approximation discussed in Section [5] (with r^j oc E-p) 

th e coeffi cient r^r is readily calculated with linearized kinetic theory for a massless 
gaPESl 

2 

rir. = 6^ , (79) 

e + V 

The kinetic theory relations A2 = —IriTT^ and A3 = are respected for the same 
reasons as the full theory. Also in the relaxation time approximation one finds, 
Ai — rjTT^. In the full kinetic theory the difference Ai — jyr^ reflects the deviation 
from the quadratic ansatz discussed Section \E\ and to a much lesser extent the 
non-linearities of the collision integral. Nevertheless the relation Ai = 7777^ almost 
holds indicating the dominance of the streaming term. Overall the relaxation time 
approximation provides a good first estimate of these coefficients in kinetic theory. 
This is important because the second order corrections to a 0-1- 1 Bjorken expansion 
vanish if Ai — tit^^ - see below. 
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From a practi cal persp ective a majority of simulations have used the full Israel- 
Stewart equationJSEnilll and treated r^r as a free parameter, varying [?y/(e + 7^)t^] 
down from the relaxation time value by a factor of two. While it is gratifying 
that higher order transport coefficients can be computed and classified, the final 
phenomenological results (see Section are insensitive to the precise value of all 
second order terms for r//s < o ^i SllSOflO l xhus, the full hydrodynamic simulations 
corroborate the estimate given in Section 14.41 for the range of validity of hydrody- 
namics. 



4.7. Summary 

We have discussed various orders in the gradient expansion of hydrodynamics. Here 
we would like to summarize these results for a + 1 dimensional Bjorken evolution. 
The equation of motion for the Bjorken expansion is 

de e + T^ 



or r 

where = t^T"^"^ is the stress te nsor at m id space-time rapidity, 77s = 0. The 
stress tensor through second order i jl3 | 25 | 3 1J 

T"'=V + (81) 

Each additional term reflects one higher order in the hydrodynamic expansion pa- 
rameter [?7/(e -|-p)t] discussed in Section We have made use of the intermediate 
results 

a'^" = diag a™, rV"") - f 0, A, ^ ^ _ 4 \ ^ ^^^^ 



3r 3r 3t 



and 



(fi^V/) - diag (^0, 



4 4 
'972'"972' 9t2 / 



(83) 



Notice in Eq. ([5T|) that there is a cancellation between the relaxation terms ^ riTT^Da 
and the visco-elastic responsJ^, Xiaa. In kinetic theory the difference TyrTr — Ai 
is determined primarily from the deviation of 5f from the quadratic ansatz (see 
Section!!]). Thus -qT-j^ is expected to be approximately equal to Ai due to the overall 
kinematics of the streaming terrrPSl, pt^Q^^f . Examining Table we see that for 
a relaxation time approximation ?7r^ = Ai and the second order corrections to a 
-fl Bjorken expansion of a conformally invariant fluid vanish! This cancellation is 
partially present for the full kinetic theory and for the strongly interacting theory. 

Fig. [14] shows the how the different orders in Eq. (|8T]) influence the evolution of 
the energy density for a conformally invariant equation of state {V — e/i = sT/A) 
and various values of 7//s. For definiteness we have used the Af = A ratios for the 
second order transport coefficients but this makes little difference since only the 
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Navier Stokes 
2nd Order 



□ 
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Fig. 14. The energy density (xt) relative to the initial energy density (xtq) for a + 1 dimensional 
Bjorken expansion. The temperature is To ~ 300 MeV and To ~ Ifm, so that 1/toTo ~ 0.66. The 
second order correction is smaller th an expected due to a cancellation between the relaxation term 
~ rjTT^Dcr and the viscoelastic terrrPD^ ^ Aicrcr. In a relaxation time approximation the second 
order correction vanishes (see text). 

combination A — rjT.^ matters - see Table 14.61 Finally we have used the estimates of 
Section [4.41 for the initial temperature To and time Tq 



Generally the effect of second order terms is small (due to the cancellation) and 
the value of the first order terms drive the correction to the ideal evolution. 

5. Kinetic Theory Description 

In Section |4] we discussed various aspects of viscous hydrodynamics as applied to 
heavy ion collisions. Since ultimately the experiments measure particles, there is a 
need to convert the hydrodynamic information into particle spectra. This section 
will provide an introduction to the matching between the kinetic and hydrodynamic 
descriptions. This will be important when comparing the hydrodynamic models to 
data in Section [6l In addition, since Section [3] discussed various calculations of 
the shear viscosity in QCD, this section we will sketch briefly how these kinetic 
calc ulations ar e performed. Good summaries of this set of steps are provided by 



In kinetic theory the spectrum of particles in a volume S is given by the Cooper- 




I^^fc l99l73ll()()l 
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Frye formulall^ 

^S^=(2iF//^^'^'^^-''-^^- ^''^ 

Note that when E is a three volume at fixed thiie, dE^ — (dV, 0,0,0), and this 
formula reduces to the traditional result. Four vectors are denoted with capitol 
letters = (i?p,p), and the equilibrium distribution function is denoted with 

-(-^-)= exp(-P.W)±l - ^''^ 

We will also use a suffix notation, Up — n{—P ■ U) and /p = f{—P ■ u). The 
distribution function obeys the Boltzmann equation 

5Jp + Vp ■ 9,/p = - / ri2^34 (/1/2 - /3/4) , (86) 
J 234 

where Wp = dEp/dp = p/Ep, and we have assumed 2^2 scattering with classical 
statistics for simplicity, fp ~ exjp(~Ep/T). The momenta are labeled as /2 — fp^ 
and /i — /pj with pi = p. The integral over the phase space is abbreviated 

r 6?p2 d^Pa <i?Pi 
234 J (27r)3 (27r)3 (2^)3 ' ^ > 

and the transition rate ri2^34 for 2 — s- 2 scattering is related to the usual Lorentz 
invariant matrix element \M\ by 

ri2-.34 = -, w ^"^J 77 r(27r)'^(5''(Pi + P2 - P3 - Pa) ■ (88) 

12-^34 (^2Ei){2E2){2E3){2Ei)^ ) \ 1^ 2 3 4) k j 

The generalization of what follows to a multi-component gas with quantum statis- 
tics is left to the references'^. 

During a viscous evolution the spectrum will be modified from its ideal form 

f = np + Sfp , (89) 

and this has important phenomenological consequences^. The modification of the 
distribution function depends on the details of the microscopic interactions. In 
a linear approximation the deviation is proportional to the strains and can be 
calculated in kinetic theory. When the most important strain is shear, the deviation 
Sf is proportional . Traditionally we parameterize the viscous correction to the 
distribution in the rest frame of the medium bj0 x(|p|) 

'5/p = -«pX(|p|)p'pV.,. (90) 

When quantum statistics are taken into account this should be written 
Sfp = -rip(l ± Tip) x(|p|)p'P^o"ij . 
where the overall minus is introduced because in the Navicr-Stokes theory tt'^" = — T^CTf" 
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Then the stress tensor in the local rest frame is 

Substituting Eq. (j90p for 6f-p and using rotational symmetry we have 

d^Tp p 
15 J (27r)3 

Thus we see that the form of the viscous correction to the distribution function 
determines the shear viscosity. 

To calculate the transport coefficients the Boltzmann equation is analyzed in 
the rest frame of a particular location x,,. In a neighborhood of this point the 
temperature and flow fields are 

u''(x,i) ~ (l,uXx,i)), T{x,t)~To + ST{x,t), (93) 

where m*(xo, t) = ST{x.o, t) = 0. The equilibrium distribution function in this neigh- 
borhood is 

n{-P ■u)^nl + n° (^^ <5T(x, t) + EI^!^^ , (94) 

where we have used the short hand notation, rt" = cxp(— i?p/To). We can now sub- 
stitute the distribution function into the Boltzmann equation and find an equation 
the Sf. The left hand side of the Boltzmann equation involves gradients, and there- 
fore only the equilibrium distribution needs to be considered. Substituting Eq. (j94p 
into the l.h.s. of Eq. ((86)) . using the ideal equations of motion 

dtu'^^- --, (95) 

dte^ -{e + V)d.,u\ (96) 



and several thermodynamic relationships 

de 

''"^dT' 



(97) 



difi/T)^—^dr + d(l] =0, (98) 



e + / T{e + V) \T 

find thalE 



dtfe + Vp ■ 9x/e = TT 



1 djW H (7,, 

3T T Tcy 2T ^ 



(99) 



The result is proportional to two strains and CTy which are ultimately respon- 
sible for the bulk and shear viscosities respectively. For a massless conformal gas 



^We have tacitly assumed that the dispersion curve E{p) does no t depend on the temperature. 
This is fine as long as we are not considering the bulk viscosity 
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we have |pp — E^, Tcv = 4e and e + p = (4/3)e. The result is that the term pro- 
portional to diU^ vanishes and consequently the bulk viscosity is zero in this limit. 
Subsequently, we will consider only the modifications due to the shear viscosity and 
refer to Section [3] for a more complete discussion of bulk viscosity. 

In Eq. ([99]) . the l.h.s. of the Boltzmann equation is evaluated at the point Xq. 
We also evaluate the r.h.s. of the Boltzmann equation at the point Xq to linear 
order in df 



o P I O O 



5f{p) , Sh 5h SU 



(100) 



/234 

In writing this equation we have made use of the detailed balance relation 

nin2Ti2^34 = n3n4ri2_»34 • (101) 

Eq. (jlOOp should be regarded as a matrix equation for the distribution function 
(5/(p). Although Sf{p) or x{p) can be determined numerically by straight for- 
ward d iscreti zation and matrix inversion, a variational method is preferred in 
practicj ^^ l ^^ i After determining 6f{p) or x(|p|) the shear viscosity can be de- 
termined from Eq. ([5^ . 

Inverting the integral equation in Eq. (|100p requires a detailed knowledge of 
the microscopic interactions. Lacking such detailed knowledge, one can resort to a 
relaxation time approximation, writing the Boltzmann equation as 

aj + .p.a./^-^-pL_i^./, (102) 

where ti^(—P- U) is a momentum dependent relaxation time. In the local rest frame 
this reduces to 

dtf + vp-d^f = ^Sf. (103) 

By fiat the correction to the distribution function is simple 



Sf = -"p^^^^'Va., . (104) 
First we consider the case where the relaxation time grows with energy 



tr{Ep) ^ Const X ^ (5/ = -Const x Up^^cjij (105) 

The form of this correction is k nown as th e quadratic ansatz and was used by all 
hydrodynamic simulations so faJSHnEIIlO]^ Substituting this form into Eq. ([92]) one 
determines that for a classical gas of arbitrary mass the constant iJ^l 

Const = ^ , with TR^Ep) cx E^ . (106) 

For a Bose or Fermi gas we have the replacement rip rtp(l ± ?ip) and Eq. (|106p 
is approximate holding at the few percent level. For a mixture of different classi- 
cal particles with one common relaxation time Eq. (|106p also holds. In practical 



May 14, 2009 21:38 WSPC/INSTRUCTION FILE hydro 



38 Derek A. Teaney 

simulations this correction is written covariantly and the phenomenological field 
TT^'' = — J/cr^ii/ is used, leading to the final result 



^/=:T7-n7^^"^X- (107) 



1 

2{e + V)T^ 

The quadratic ansatz may seem arbitrary, but it is often a good model of coUi- 
sional energy loss and weak scattering. For instance an analysis of the leading-log 
Boltzmann equation (along the lines of Eq. (I100[) ) shows that the quadratic ansatz 
describes the full results to 10-15% accuracj^^. However, this agreement is in part 
an artifact of the leading-log, or soft scattering, approximation. For example, in 
th e lea ding-log plasma the energy loss of a "high" pT quark from the bath is given 



bjUnsi 

From this formula we see that the energy loss is constant at high momentum and 
therefore the relaxation time scales as tji (x p in agreement with Eq. (|105p . In 
reality Eq. (jlOSp is decidedly wrong at large momentum where radiative energy 
loss becomes increasingly significant and can shorten the relaxation time. Indeed 
when coUinear radiation is included in the Boltzmann equation the quadratic ansatz 
becomes increasingly poor Section[S]we will consider a relaxation time which 

is independent of energy as an extreme alternative 

TR oc Const , (109) 

and explore the phenomenological consequences of this ansatz. 

As discussed in Section [51 the differential elliptic flow V2{pt) is sensitive to the 
form of these corrections, while the integrated V2 is constrained by the underly- 
ing hydrodynamic variables, and is largely independent of these details. This last 
remark should be regarded with caution as it has not been fully quantified. 



6. Viscous Hydrodynamic Models of Heavy Ion Collisions 

At this point we are in a position to discuss several viscous hydrodynamic models 
which have been used to confront the elliptic flow data. To initiate discussion, 
we show simulation results for V2{pt) from Luzum and Romatschke in Fig. 1151 
Comparing the simulation to the "non-flow corrected" data for pT ^ 1.5 GeV, we 
can estimate an allowed range for the shear viscosity, ry/sR:! 0.08 ^0.16. Below we 
will place this estimate in context by culling figures from related works. 

A generic implementation of viscous hydrodynamics consists of several parts: 

(1) At an initial time Tq the energy density and flow velocities are specified. For 
Glauber initial conditions, one takes for example 

d^' 
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Glauber 
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Fig. 15. Figure from Ref.ED which shows how elliptic flow depends on shear viscosity. The theory 
curves are most dependable for px ^ 1.5 GeV and should be compared to the "non-flow corrected" 
data. The Glauber and CGG initial conditions have different eccentricities as described in the text. 



where the overall constant is adjusted to reproduce the multiplicity in the event. 
The simulations assume Bjorken boost invariance with the ansatz 

e(r,x^,77) = e(r,x^) , (111) 

U^{t,^^,7^) = {U\u\uy,u'^) = (u"(T,Xi), U^{t,X^), Uy{T,^^), 0) . (112) 

In cartesian coordinates = u'^sinh(77s) and m* = u"^ cosh(7ys). The calcula- 
tions typically assume zero transverse flow velocity at the initial time Tq 

M"(ro,Xi) =u'^(ro,Xi) = 0, u"(To,Xi) = l. (113) 
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The strains are taken from the Navier stokes theory for example 

^''''(r„,x^) -diag(7r-^7^-^^^^r27r'"') = (^0, -j^^ , (114) 

and reflect the traceless character of shear stress. 

(2) The equations of motion are then solved. Viscosity modifies the hydrodynamic 
variables, T and u^, and also modifies off diagonal components of the stress 
tensor through the viscous corrections tt^". 

(3) A "freezeout" condition is specified either by specifying a freezeout temper- 
ature or a kinetic condition. During the time evolution a freezeout surface is 
constructed. For instance the freezeout surface in Fig.[Tn]is the space-time three 
volume S where Tfo ~ 150 MeV. 

(4) Finally, in order to compare to the data, particle spectra are computed by 
matching the hydrodynamic theory onto kinetic theory. Specifically, on the 
freezeout surface final particle spectra are computed using Eq. (|125|) . Roughly 
speaking this "freezeout" procedure is equivalent to running the hydro up to 
a particular proper time Tf or temperature Tf and declaring that the thermal 
spectrum of particles at that moment is the measured particle spectrum. 

There are many issues associated with each of these items. The next subsections 
will discuss them one by one. 



6.1. Initial Conditions 

First we note that the hydrodynamic fields are initialized at a time tq — 1 fm/c, 
which is arbitrary to a certain extent. Fortunately, both in kinetic the ory and 
hydrodynamics the final results are not particularly sensitive this value 1 ^4 1 31 1 
Also, all of the current simulations have assumed Bjorken boost invariance. While 
this assumption should be relaxed, past experience with ideal hydrodynamics shows 
that the mid-rapidity elliptic flow is not substantially modified . Above we have 
discussed one possible initialization of the hydro which makes the energy density 
proportional to the number of binary collisions, e.g. the Glauber curves of Fig. 1151 
Another rea son able option is to make the entropy proportional to the number of 
participants 

s(to,x_l) oc . (115) 

da; ay 

As a limit one can take the CGC model discussed in Section[51 Finally it is generally 
assumed that the initial transverse flow is zero 

Ua:(To,X^) = Wj^(to,X^) = 0. (116) 

This assumption should probably lifted i n fut ure calculations and a more reasonable 
(but still small) estimate is given in 

Examining Fig. [12] and Fig. [7] we see see that there is a significant and pre- 
dictable linear dependence on the eccentricity. When extracting the shear viscosity 
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from the data, this uncertainty in the eccentricity leads to a factor of two uncer- 
tainty in the final results for the shear viscosity. As emphasized in Section [2l the 
CGC model should be thought of as an upper limit to the anisotropy that can 
be produced in the initial state. Therefore, the uncertainty in -q/s is probably not 
as large as dispersion in the curves would indicate. In ideal hydrodynamics, the 
spread in V2{pt) resultin g fro m the different initializations specified by Eq. (jllSp 
and Eq. (|110p was studiedly, and is small compared to the difference between the 
Glauber and CGC curves in Fig. [151 

Once the initial conditions for the temperature and the flow velocities are speci- 
fied, the off diagonal components of the stress tensor tt^^'^ are determined by the spa- 
tial gradients in T and . To second order this is given by Eq. ([GS]) and there is no 
ambiguity in this result. Time derivatives may be replaced with spatial derivatives 
using lower order equations of motion to second order accuracy. In the phenomeno- 
logical theory tt^" is promoted to a dynamical variable. Clearly the appropriate 
initial condition for this variable is something which deviates from —rja^'^ by sec- 
ond order terms. However, the extreme choice tt^'^ = was studied to estimate how 
initial non-equilibrium effects could alter the final results. This is just an estimate 
since the relaxation of these fields far from equilibrium is not well captured by 
hydrodynamics. On the other hand, comparisons with full kinetic theory simula- 
tions show that the Israel-Stewart model does surprisingly well at reproducing the 
relaxation dynamics of the full simulatioiJ^. From a practical perspective, even 
with this extreme choice tt^" = 0, the stress tensor relaxes to the expected form 
tt'"' = —rjc^'^ relatively quickly. The result is that V2{pt) is insensitive to the 
different initializations of tt^"^. This can be seen clearly from Fig. [TBI 

6.2. Corrections to the Hydrodynamic Flow 

Once the initial conditions are specified, the equations of motion can be solved. 
First we address the size of the viscous corrections to the temperature and flow 
velocities. The magnitude of the viscous corrections depends on the size of the 
system and the shear viscosity. Fig. [T7l shows a typical result for the AuAu system. 
As seen from the figure the effect of viscosity on the temperature and flow velocities 
is relatively small, of order ^ 15% for 77/s ~ 0.2. 



An explanation for this result is the followin g^ 'I -. In the first moments of the 
collision the system is expanding in the longitudinal direction and the pressure in 
the longitudinal directions is reduced 



At first sight, this means that the system cools more slowly and indeed this is 
initially true. However, since the shear tensor is traceless there is an increase in the 
transverse pressure which is uniform in all directions 




^Sjiijj) _ p _ _^ 



(117) 



T- 



6 T 




(118) 
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Fig. 16. Figure from studying the independence of the final results on the initialization of 

tt'^". Note that due to different metric and symmotrzation conventions the Navier Stokes limit is 
2r]a"^" rather than —rja'^'^ adopted here. 




Fig. 17. A central AuAu simulation with an ideal gas equation of state p = e/3 and rj/s = 0.2 
which compares the viscous and Euler evolutionl07 . The left fi gure shows the energy density (xr) 
for different times. The right figure shows the velocity for different times. 



and which ultimately increases the radial flow. Since the radial flow is larger in the 
viscous case, the system ultimately cools faster. 

Having discussed the dependence of T and u'^, we turn to a qu antity which 
largely dictates the final elliptic flow - the momentum anisotropy"^. The momen- 
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Fig. 18. Figure from Ref.ESl comparing the development of the flow anisotropy (Eq. I I119I I) in 
viscous hydrodynamics relative to the ideal hydrodynamics. The lower band of curves are all rep- 
resentative of viscous hydro and differ only in how the second order corrections are implemented. 
The anisotropy differs from ideal hydro because the anisotropy involves the viscous difference, 



turn anisotropy is defined a^ 

, _ (T^^ - Tyy) _ J (T^^ - Tyy) 

Fig. 1 181 illustrates how this momentum anisotropy increases as a function of time in 
the CuCu and AuAu systems. Although the flow fields T(xj_, r) and m^(xj_, r) are 
quite similar between the ideal and viscous cases, the ideal anisotropy is significantly 
reduced by viscous effects. The reason for this reduction is because the viscous stress 
tensor anisotropy, T^^ — Tyy, involves the difference 



lixx _ ^yy ^ 



in addition to the temperature and flow velocities. This additional term is ultimately 
responsible for the deviation of between the ideal and viscous hydrodynamic 
calculations. At later times there is some modification of e'p, due to the flow itself, 
but this is dependent on freezeout. The deviation AH — H^^ — ipy will have 
important phenomenological consequences in determining the viscous correction to 
the elliptic flow spectrum. 



' Note there is a misprint in the original definition of in Ref.lSU. Eq. (3.2) of that work used 
a double bracket notation indicating an energy density weight which should only be a single 
bracke t as above. This double bracket definition (rather than Eq. II119II ) was used subsequently 
in RefPSI which is why the associated curve differs from other recent works | 31 | 30 | Pm-^her the 
definition should perhaps include a factor of so that the integrals would be boost invariant, i.e. 
f dS^u'' (T^^ - Tyy) = tA-Qs f d^xu°(T^== - Tyy). However, typically the integral is dominated 
near the center where ~ 1. An alternative definition Sf^^^^^ is suggested in Section 16.51 
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6.3. Convergence of the Gradient Expansion 

Fig. [18] also compares the simplified Israel-Stewart equation Eq. ([74]) to the full 
Israel-Stewart equation Eq. (|75p as a function of the relaxation time parameter r^. 
The result supports much of the discussion given in Section 14.61 In the roughest 
approximation neither the simplified Israel-Stewart equation nor the full Israel- 
Stewart equation depend on the relaxation time parameter t^. When an ideal gas 
equation of state is used the dependence on is stronger especially for the sim- 



plified Israel-Stewart equatiorPQl. Note also that the dependence on is stronger 
in the smaller CuCu system than in AuAu. However, in the conformally invariant 
full Israel-Stewart equation the dependence on is negligible, indicating that the 
result is quite close to the first order Navier Stokes theory. The more rapid con- 
vergence of the gradient expansion in conformally invariant fluids is a result of the 
fact that the derivatives in the conformal case come together as 

We have selected one figure out of manjEnHOEI], xhe result of the analysis is that 
the fiow fields and V2{pt) are largely independent of the details of the second order 
terms at least for 77/5 < 0.3. For this range of parameters, hydrodynamics at RHIC 
is an internally consistent theory. 

6.4. Kinetic Theory and Hydrodynamic Simulations 

Clearly viscous hydrodynamics is an approximation which is not valid at early 
times and near the edge of the nucleus. This failure afflicts the current viscous 
calculations at a practical level right at the moment of initialization. For instance, 
the longitudinal pressure 

T'-' EET^T'^" = V --— , (120) 
3 To 

eventually becomes negative near the edge of the nucleus indicating the need to 
transition to a kinetic descriptiorff^. (Note that V <x while cx , so at 
sufflciently low temperatures the viscous term is always dominant regardless of the 
magnitude of jy/s.) The current calculations simply limit the size of this correction 
through the phenomenological Israel-Stewart model. For example, one approach 
would be to take 



T^T"" = 7^ + T^n"" , (121) 



with 



,.n- J"^^ while 4/3,/ro< 0.9 _ ^^^^^ 
1—0.97' otherwise 

This ad-hoc flx is clearly not nice and points to the larger problem of freezeout 
which is difficult to address with hydrodynamics itself. 
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Freezeout is the colloquial term for the transition from a hydrodynamic to a 
kinetic regime and is impossible to separate cleanly from the viscosity itself in a 
realistic nucleus-nucleus collision. Clearly as the shear viscosity is made smaller and 
smaller, a larger and larger space time volume is described by hydrodynamics. To 
estimate the size of the relevant space time region we remark that hydrodynamics 
is valid when the relaxation time is much smaller than the inverse expansion 
rate, t/jC^^m'^ <^ 1. Therefore, in the simulations one can estimate the region of 
validity by monitoring the expansion rate relative to the relaxation timJi^SIIini 
Specifically, freezeout is signaled when 



^d,u^^ ^ \ . (123) 



-d^u^' ^ TRd^uf' ^ - . (124) 



This combination of parameters can be motivated from the kinetic theory estimates 
KD. The pressure isV^e (v^^^), with (f^h) typical quasi-particle velocity. The 
viscosity is of order 77 ^ e (wth) '''R with th the relaxation time. Thus hydrodynamics 
breaks down when 

Fig. [in] estimates the space-time region described by viscous hydrodynamics. Ex- 
amining this figure we see that the time duration of the hydrodynamic regime is a 
relatively strong function of 77/s at least for a conformal gas. In reality the behavior 
of the shear viscosity near the transition region will control when the hydrodynam- 
ics will end. 

Clearly the surface to volume ratio in Fig.[T9]is not very small. Hydrodynamics is 
a terrible approximation near the edge. There is a need for a model which smoothly 
transitions from the hydrodynamic regime in the center to a kinetic or free streaming 
regime at the edge. Near the phase transition kinetic theory may not be a good 
model for QCD, but it has the virtue that it gracefully implements this hydro to 
kinetic transition. Although the interactions and the quasi-particle picture of kinetic 
theory may be decidedly incorrect, this is unimportant in the hydrodynamic regime. 
In the hydrodynamic regime the only properties that determine the evolution of the 
system are the equation of state, V{e), and the shear viscosity and bulk viscosities, 
77(e) and C{^)- III the sense that kinetic theory provides a reasonable guess as to how 
the surface to volume ratio influences the forward evolution, these models can be 
used to estimate the shear viscosity and the estimate may be more reliable than the 
hydrodynamic models. A priori one should demand that the kinetic models have the 
same equation of state and the same shear viscosity as expected from QCD, 77 cx 
For instance in a kinetic m odel of m assless particles with a constant cross section 
ctq (such as studied in the shear viscosity scales as 77 oc T/ctq. This 

difference with QCD should be kept in mind when extracting conclusions about 
the heavy ion reaction. Further, many transport models conserve particle number, 
which is an additional conservation law not inherent to QCD; this also changes the 
dynamics. Keeping these reservations in mind we examine Fig.[20]from Ref.QJJ. This 
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Fig. 19. This figure estimates for a conformal gas with equation of s tate p = e/3 and con- 
stant rj/s the space time region described by viscous hydrodynamics The contours are where 
(jj/'P) dfiU'^ = 0.6 for different values of rj/s. For the smallest value rj/s = 0.05 the system freezes 
out at a time of order ~ 40 fm. This unrealistically long time reflects the conformal nature of the 
gas as discussed in Section 14.51 For comparison we have shown the (rj/'P) d/j.u'^ = 0.225 contour 
for r//s = 0.05. 

figure shows a promising comparison between kinetic theory with a constant cross 
section (ctq = Const) and a viscous hydrodynamic calculation with rj cx T/ao- The 
case with a cx t^/'^ will not be discussed in this review, but is an attempt to mimic 
a fluid which has rj oc T^. 

What is exciting about this figure is the fact that the hydrodynamic conclu- 
sions are largely supported by the results of a similar kinetic theory. This gives 
considerable confidence that surface to volume effects are small enough that the 
hydrodynamic conclusions presented in Fig. [15] are largely unchanged by particles 
escaping from the central region. More formally, the opacity is large enough to 
support hydrodynamics. 

There have been other kir ietic calcu lations which are working towa rds extracting 
ry/s from the heavy ion data I12L 32 | 33 | jj^ particular Refs.^^^S^^H^IIISI ^ged a kinetic 
theor y irn plementation of 2 — > 2 and 2 — > 3 interactions motivated by weak coupling 
QCD ISSI^ The simulation also calculates the Debye scale self consistently, i.e. in 
equilibrium one sets mj^ cx g^T^ and T changes with time. Out of equilibrium 
this mass scale is determined from the momentum distribution of particles. 
Consequently this model respects the symmetry properties of high temperature 
QCD, i.e. the model has rj (x and does not conserve particle number. For the 
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Fig. 20. Elliptic flow from Ref.UDfor a massless classical gas with a constant cross section in kinetic 
theory (Transp. a = const) and viscous hydrodynamics (Hydro cr = const). The a oc r^/^ case is 
not discussed in this review but is an effort to simulate a gas with rj oc T^. 



model parameter — 0.3 <-!■ 0.5 (which is only schematically r elated to the running 
coupling) the shear to entropy ratio is r]/s = 0.16 <-> 0.08 | H3 | 112| rpj^^ model 
(known as BAMPS) is conformal and never freezes out as discussed in Section l4?5l 
The current implementation of BAMPS simply stops the kinetic evolution when the 
energy density reaches a critical value, Cc — 0.6 ^ l.OGeV/fm'^. This is an abrupt 
way to introduce a needed scale into the problem and schematically approximates 
the rapid variation of the shear viscosity in this energy density range. Fig. l21l shows 
the time development of elliptic flow in this model which can reproduce the observed 
flow only for rj/s = 0.16 ^ 0.08. The time development of V2 seen in Fig. [21] shows 
that it is very difhcult to separate precisely the shear viscosity in the initial stage 
from the freezeout process controlled by Cc. 

Clearly the transition from a hydrodynamic regime to a kinetic regime is im- 
portant to clarify in the future. In the meantim e most hy drodynamic groups have 
invoked an ad-hoc freezeout prescription. In Refs.!^ ^ ! *^^ ! ^ ! the hydrodynamic codes 
were run until a typical freezeout temperature, Tfo — 150 MeV. Technically the 
freezeout surface is constructed by marching forward in tim e and triangulating 
the space-time surface with constant temperature. In a surface of constant 

particle density was chosen n ~ 0.365/fm^, and in Rem the chosen surface was 
motivated by the kinetic kinetic condition in Eq. (|123p . Ideally this could be im- 
proved by dynamically coupling the hydrodynamic evolution to a kinetic description 



I I I I I I I I I I I I I I I I I I I I I I I _ 

■ Transp. a = const. — 
n Transp. .CTocT^/^i^ ^ 




May 14, 2009 21:38 WSPC/INSTRUCTION FILE hydro 



48 Derek A. Teaney 



0.09 r- 

0.08 



^2 0.07 - 




part 



Fig. 21. (Top) The development of elliptic flow V2 as a function of time in the BAMPS modeH 
(Bottom) The final elliptic flow as a function of centrality. The shear viscosity t o entropy ratio 
rj/s corresponding to the model parameter as = 0.3 <-> 0.6 was estimated in and is 

r]/s = 0.16 <-> 0.08. The evolution is stopped when the energy density reaches a critical value of 
ec- 



or by simulating the entire event with a kinetic model which closely realizes the 
equation of state and transport coefficients used in the hydrodynamic simulations. 



6.5. Particle Spectra 

Finally we turn to the particle spectra in viscous hydrodynamics. Ideally the sys- 
tem would evolve through the approximate phase transition down to sufficiently 
low temperatures where the dynamics could be described either with viscous hy- 
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drodynamics or with the kinetic theory of a Hadron Resonance Gas (HRG). In 
reahty this does not seem particularly likely since the system is already expanding 
three dimensionally and the scales are approximately fixed (see Section |4?5|) . The 
estimates of the shear viscosity to entropy ratio in a hadronic gas are reliable for 
T <, 130 MeV and do not support this optimistic picture (see Section [3]). In seems 
quite unlikely that there is equilibrium evolution in the HRG below a temperature of 
T ~ 150 MeV. Clearly the dynamics is extremely complex during the quark-hadron 
transition which takes place for an energy density of e ~ 0.5 ^ 1.2GeV/fm^. In 
this range, the temperature changes by only AT ~ 20 MeV. However, the hydro- 
dynamic simulations evolve this complicated region for a significant period of time, 
T ~ 4 fm <-!■ 7 fm. This transition region can be seen from the inflection in the AuAu 
plots in Fig. [HI 

The pragmatic approach to this complexity is to compute the quasi-particle 
spectrum of hadrons at a temperature of T ~ 150 MeV. Since the HRG describes the 
QCD thermodynamics well, this pragmatism is fairly well motivated. The approach 
conserves energy and momentum and when viscous corrections are included also 
matches the strains across the transition. In ideal hydrodynamics simulations the 
subsequent evolution of the hadrons has been followed with hadronic cascade models 
[22 19 18^ rpj^g result of these hybrid models is that the hadronic rescattering is 
essentially unimportant for the V2{pt) observables presented here. 

Technically, the procedure is the following: along the freezeout surface the spec- 
trum of particles is computed with the Cooper-Frye formula 

where a labels the particle species, the distribution function is, 

n-P ■ u) = n^i-P ■ u/T) + 5r{-P ■ u/T) , (126) 

and da labels the spin-isospin degeneracy factor for each particle included ( see 
Section [5]). In practice, the Boltzmann approximation is often sufficient. In Ref.EIl 
all particles were included up to mass of mros < 2.0 GeV and then subsequently 
decayed. In other works a simple single species gas was used t o study various aspects 
of viscous hydrodynamics divorced from this complex realitjff^^lil 

All of the viscous models used the quadratic ansatz discussed in Section [5l 
writing the change to the distribution function of the a-th particle type as 

r^n'^ + Sr , (127) 

with (5/" given by 

Sr - ^—r^n'^n ± n'')PfP''Tr,,^ . (128) 

Before continuing we review the elements that go into a complete hydrodynamic 
calculation. First initial conditions are specified (see Section WT^ ; then the equa- 
tions are solved with the viscous term (see Section [6T2|) ; after this a freezeout surface 



E- 
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Fig. 22. V2(pt) based on Refs.HSHSI showing the V2{pt) with (fo + Sf) and and without (/o only) 
the viscous modification of the distribution function. The result depends to a certain extent on 
freezeout and the freezeout temperature here is Tfo = 130 MeV. 



is specified (see Section for the hmitations of this); finally we compute spectra 
using Eq. (|125p and Eq. p28p . This particle spectra can ultimately be compared 
to the observed elliptic flow. With this oversight we take a more nuanced look at 
Fig. [151 

To separate the viscous modifications of the flow variables {T and u^) and the 
viscous modifications of the distribution function, we turn to Fig. 1221 Examining 
this figure we see that a significant part of the corrections due to the shear viscosity 
are from the distribution function rather than the flow. Although the magnitude 
of the flow modifications depends on the details of freezeout, this dependence on 
(5/ is the typical and somewhat distressing result. We emphasize however that it 
is inconsistent to drop the modifications due to Sf. The result of dropping the Sf 
means that the energy and momentum of the local fluid cell T^^'^Uu is matched 
by the particle content, but the off diagonal strains tt^'^ are not reproduced. The 
assumption underlying the comparison of viscous hydrodynamics to data is that 
the form of these off diagonal strains is largely unmodified from the Navier Stokes 
limit during the freezeout process. As the particles are "freezing out" this is a 
reasonable assumption. However, since these strains are only partially constrained 
by conservation laws, this assumption needs to be tested against kinetic codes as 
already emphasized above. This freezeout problem is clearly an obstacle to a reliable 
extraction of ry/s from the data. 

Although the dependence on 5 J in V2{pt) is undesirable, the viscous modifi- 
cations of the integrated elliptic flow V2 largely reflects the modifications to the 
stress tensor itself. The observation is that the stress anisotropy e' (see Eq. (|119p ) 
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Fig. 23. Dependence of elliptic flow versus some centrality from The lines show the results 

of viscous hydrodynamics and the colored squares show the anisotropy of the stress tensor 0.5 Cp 
(see Eq. II119I I) for different values of ri/s. The sensitivity to the quadratic ansatz is estimated in 
the text and corresponds to half the difference between the red (r)/s = 0.08) and blue {rj/s = 0.16) 
curves. The heuristic rule V2 ~ 0.5 e' is motivated in the text. 



determines the average flow according to a simple rule of 

thumlJSIEI] 

V2^\e',. (129) 

Fig. [221 shows e'p and V2 as a function of centrahty. The figure supports this rule 
and suggests that sufficiently integrated predictions from hydrodynamics do not 
depend on the detailed form of the viscous distribution. 

To co rro borate this conclusion we turn to an analysis or iginally presented by 
OllitraulJISl and subsequently generalized to the viscous cas e ^ I First we pa- 
rameterize the single particle spectrum dN/dpT with an exponential and V2{pt) as 
linearly rising, i.e. 

— 4^~Ce-P^/^, and V2 (x pr ■ (130) 

PT dpT 

With this form one finds quite generally that the p^ weighted elliptic flow is twice 
the average V2 

The p^ weighted elliptic flow has a much closer relationship to the underlying 
hydrodynamic variables. Indeed we will show how this simple rule of thumb arises 
and that it is largely independent of the details of 5f. 
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To this end, we evaluate the sphericity tensor which will have a simple relation- 
ship to A2 

gi^up ^ g^up ^ g^up ^ j _^P_p^^p,^pp [„(„p . ^) ^ gf^_p . J,)] (132) 

Then using the Cooper-Frye formula (with dS^ — (rdf^s d^x^, 0, 0, 0)) the asym- 
metry at any given moment in proper time is 

The sphericity tensor consists of an ideal piece and a viscous piece 

f->n + 6f, S'^'P -> ST" + S'^.''" . (134) 

First we consider the ideal piece and work in a classical massless gas approximation 
for ultimate simplicity. The tensor is a third rank symmetric tensor and can be 
decomposed as 

gpi^P ^ + B{T) {u^'g^P + perms) . (135) 

Here A{T) and B{T) are thermodynamic functions and are given by 

A f d^V , . , . 

(2^y"^-(^+^)^- ^'^'^ 

For Bose and Fermi gases this relation between A{T), B{T), and the enthalpy is 
approximate. Thus the ideal piece of the sphericity tensor is largely constrained 
by thermodynamic functions. The viscous piece is largely constrained by the shear 
viscosity. As discussed in Section [5] we parameterize the Sf with x(p) 



x{p) 
{p-uy 

While the precise form of the viscous correction x{p) depends on the details of the 
microscopic interactions, it is constrained by the shear viscosity 

Substituting the viscous parameterization into the definition of the sphericity we 
find 

goxx ^ -C{T) [itV^^ + 2m^ct^°] , (139) 
~ -C(r)MV"" , (140) 



where C(T) is 
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For simplicity we have assumed that the flow is somewhat non-relativistic so that 
yXj^xa O^v^) compared to u^tt^^. To get a feeling for how sensitive the results 
are to the quadratic ansatz we will work with a definite functional form 

X{p) = Const X p^-" . (142) 

In a relaxation time approximation discussed in Section [5l a — Q corresponds to 
a relaxation time which increases with linearly p while a = 1 corresponds to a 
relaxation time independent of p. Substituting this ansatz we find 

C(T) ~ (6 - a)T?7 , where < a< 1 . (143) 

Having assembled the ingredients we can write down an approximate formula 
for A2 

^ J d^xj^ Tu° [(e + P)(m^m^ - uyuy) + (1 - f ) (tt^"' - tt?'?')] 

^ ~ / d2x^ TwO [(e + 'P){u''u^ + uVuV) + (e + 'P)/'i+{l- f ) (vr^^ + tt^)] ' 

(144) 

with < a < 1. This is the desired formula which expresses the observed elliptic 
flow in terms of the hydrodynamic variables. To reiterate, the coefficient a changes 
the functional form the viscous distribution function and < a < 1 is a reasonable 
range — a = is the usual quadratic ansatz. It is useful to compare this formula 
to the definition of 

/d^xx [(e + P)(u"^u"=-u2'u?') + (7r^^-7r?'2')] 
~ /d2x_L [(e + 'P)(u^u^ + M^'u!') + 27^ + (tt^^ + TT?'!^)] ■ ^ ' 

Thus while is not exactly equal to the A2 of Eq. p44p . it is close enough to 
explain the heuristic rule, 2v2 — e!^. 

The overall symmetries and dimensions of the sphericity tensor suggests a defi- 
nition for an analogous quantity in hydrodynamics 



QilVp _ rp 

■^hydro " ^ 



(e + V)u^'v"u'' + ^(e + V){u^g'''' + perms) + (u'^tt'^" + perms) 



(146) 

In fact in a classical massless gas approximation with the quadratic ansatz, the 
analysis sketched above shows that 

The important point for this review is that from Eq. (|144p we see that the integrated 
elliptic flow is relatively insensitive to the quadratic form for the viscous distribution 
function. More specifically the uncertainty is ^ 15% of A^^'^'^ — ^^'scous^ j g_ about 
half the difference between the blue (77/5 = 0.16) and red curves (77/3 = 0.08) 
in Fig. [23l However, the estimate shows that changing the quadratic ansatz will 
move the curves systematically higher increasing the preferred value of 7]/s to a 
certain extent. The quadratic ansatz (which has been used by all hydrodynamic 
simulations) can be a poor approxi mation when coUinear emission processes are 
included in the Boltzmann descriptiorff^. 
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Clearly addressing more completely the uncertainties associated with the parti- 
cle content and microscopic interactions in Fig.[23]is a task for the future. Neverthe- 
less, it does seem that sufficiently integrated quantities will reflect rather directly 
the bulk properties of the hydrodynamic motion in a way that can be quantified. 

7. Summary and Outlook 

The elliptic flow data presented in Fig. [5] and Fig. [S] provide strong evidence for 
hydrodynamic evolution in a deconfined phase of QCD. It is difficult to think of 
other QCD based mechanisms which could explain the measured elliptic flow. This 
is because the nuclear geometry has a size Rau ~ 4fm, which is large compared 
to the momentum scales in the problem p ~ 0.6 GeV. Given this large size, the 
response of the nuclei to this initial geometry must develop over a relatively long 
time, T ~ 4fm/c. Hydrodynamics is the appropriate framework to describe the 
collective motion over these long time scales. 

We have discussed three advances which have corroborated the hydrodynamic 
interpretation of the flow. The first advance is experimental and a brief overview 
is provided in Section [21 These measurements now show quite clearly that the 
hydrodynamic response is from a deconfined phase (see Fig. [5]). Furthermore, the 
measurements also show that the flow decreases in smaller systems in a systematic 
way (see Fig. [S|). 

The second advance is in viscous hydrodynamics. There has been important 
conceptual progress in understanding hydrodynamics beyond the Navier-Stokes or- 
der (see Section 14.61) . The important result of this analysis (see Fig. [14)) is that for 
a Bjorken expansion there are generic kinematic cancellations between the second 
order terms which reduce their relative importance. In fact in a -I- ID Bjorken 
expansion of a conformal gas in the relaxation time approximation, the second 
order corrections vanish, while in other more general kinetic theories the second or- 
der corrections almost vanish. This understanding of second order hydrodynamics 
has spurred additional progress in viscous simulations of the heavy ion event (see 
Section [6]). It is satisfying that viscous hydrodynamics simulations are in better 
agreement with the data than ideal hydrodynamics and naturally explain several 
trends in the elliptic flow. For example, viscosity explains the fall off of V2 at high 
momentum (see Fig. [22l and Fig. [5]) and the fall off of V2 in more peripheral colli- 
sions (see Fig. [Hand Fig. [23]). Many of these trends were previously reproduced by 
transport modelJ^] and it is exciting to see that they are now reproduced with a 
macroscopic approach. 

The final advance has been in kinetic theory simulations. Kinetic theory simu- 
lations smoothly interpolate between the hydrodynamic regime and free streaming. 
Since hydrodynamics is universal {i.e. it only depends on V{e) and ?y(e)), the kinetic 
theory results can be used to estimate how still higher gradients would influence 
the hydrodynamic conclusions. Ideally the results of these simulations would be 
largely independent of the details of the interactions. The kinetic simulations have 
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demonstrated the abihty to reproduce the viscous hydrodynamics (see Fig. 1^01) ). 
validating the hydrodynamic assumption that the dynamics near the edge does not 
significantly influence the time development of elhptic flow. In addition, various ad- 
vances have allowed these kinetic codes to simulate a fluid which does not conserve 
particle number and which has a shear viscosity, r] oc (see Fig. [2T|) . Thus, the 
estimates of rj/s from these simulations are complementary to the macroscopic ap- 
proach, and must be accepted even if theoretical prejudice rejects the microscopic 
details. 

Examining the figures in this review, we see t hat there is consensu s between 
hydrodynamic and kinetic models on several point j ^ l ^ l ^ l ^ l ^ l ^ l -*^*^ ! ^^ ! -*^^ 

• Shear viscosity is needed to reproduce trends in the data. This consensus has 
yet to translate to an agreed upon lower bound on rj/s but probably will in the 
not too distant future. 

• It is impossible to reproduce the elliptic flow if the shear is too large: 

?7/s>0.4. (148) 

This bound is quite safe and has been found by all groups which have tried to 
reproduce the observed flow. 

• The preferred value of rj/s is (see Fig. [TSl and Fig. [2T|) 

v/s ~ (1 ^ 3) X i- . 
To reduce these constraints further several items need to studied and quantifled. 

• At the end of the collision kinetic assumptions need to be made. The uncer- 
tainties involved with the particle content and the quadratic ansatz to the 
distribution need to be quantified to a much greater extent than has been done 
so far. The theoretical motivations and limitations of the quadratic ansatz have 
been discussed in Section [5l In Section we have taken the nascent steps to 
quantify these uncertainties with an independent analysis of Fig. [23] due to 
Luzum and Romatschke. 

• The nucleus-nucleus collision terminates as the system enters the transition 
region and starts expanding three dimensionally. Some estimates for this tran- 
sition process have been given in Section [4751 Clearly the transition region sets 
a very definite scale e ~ 0.6 ^ 1.2GeV/fm'^ which can not be ignored. It is 
difficult to separate the rapid 77(e) dependence in this region from the shear 
viscosity well into the QGP phase (See Fig. [2T|) . This scale influences the de- 
pendence of V2 on centrality and hampers an extraction of rj/s from this size 
dependence. It will be important to categorize in a model independent way {e.g. 
77(e)) how this scale influences the flnal flow and the flnal estimates of 77/s. 

Finally, while this review has focused squarely on elliptic flow, the short trans- 
port time scales es tim ated here have implication s for a large n umb er of othe r obsery- 
ables - energy losJ^^, the ridge and Mach coneJU^m^'II^, heavy quarks H^^^, 
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and many more. Ultimately these observables will provide a more complete picture 
of strongly coupled dynamics near the QCD phase transition. 
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